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ABSTRACT 

In this series of lectures, aimed for non-specialists, I review the consider- 
able progress that has been made in the past decade in understanding how 
galaxies form and evolve. Complementing the presentations of my theoretical 
colleagues, I focus primarily on the impressive achievements of observational 
astronomers. A credible framework, the ylCDM model, now exists for inter- 
preting these observations: this is a universe with dominant dark energy whose 
structure grows slowly from the gravitational clumping of dark matter halos 
in which baryonic gas cools and forms stars. The standard model fares well 
in matching the detailed properties of local galaxies, and is addressing the 
growing body of detailed multi- wavelength data at high redshift. Both the 
star formation history and the assembly of stellar mass can now be empir- 
ically traced from redshifts z ~6 to the present day, but how the various 
distant populations relate to one another and precisely how stellar assembly 
is regulated by feedback and environmental processes remains unclear. In the 
latter part of my lectures, I discuss how these studies are being extended to lo- 
cate and characterize the earliest sources beyond z ~6. Did early star- forming 
galaxies contribute significantly to the reionization process and over what pe- 
riod did this occur? Neither theory nor observations are well-developed in this 
frontier topic but the first results are exciting and provide important guidance 
on how we might use more powerful future facilities to fill in the details. 
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1 Role of Observations in Cosmology Galaxy 
Formation 

1.1 The Observational Renaissance 

These are exciting times in the field of cosmology and galaxy formation! To 
justify this claim it is useful to review the dramatic progress made in the 
subject over the past ~25 years. I remember vividly the first distant galaxy 
conference I attended: the lAU Symposium 92 Objects of High Redshift, held in 
Los Angeles in 1979. Although the motivation was strong and many observers 
were pushing their 4 meter telescopes to new limits, most imaging detectors 
were still photographic plates with efficiencies of a few percent and there was 
no significant population of sources beyond a redshift of z=0.5, other than 
some radio galaxies to z ~1 and more distant quasars. 

In fact, the present landscape in the subject would have been barely rec- 
ognizable even in 1990. In the cosmological arena, convincing angular fluc- 
tuations had not yet been detected in the cosmic microwave background nor 
was there any consensus on the total energy density Htot- Although the 
role of dark matter in galaxy formation was fairly well appreciated, neither 
its amount nor its power spectrum were particularly well-constrained. The 
presence of dark energy had not been uncovered and controversy still reigned 
over one of the most basic parameters of the Universe: the current expansion 
rate as measured by Hubblc's constant. In galaxy formation, although evolu- 
tion was frequently claimed in the counts and colors of galaxies, the physical 
interpretation was confused. In particular, there was little synergy between 
observations of faint galaxies and models of structure formation. 

In the present series of lectures, aimed for non-specialists, I hope to show 
that we stand at a truly remarkable time in the history of our subject, largely 
(but clearly not exclusively) by virtue of a growth in observational capabili- 
ties. By the standards of all but the most accurate laboratory physicist, we 
have 'precise' measures of the form and energy content of our Universe and a 
detailed physical understanding of how structures grow and evolve. We have 
successfully charted and studied the distribution and properties of hundreds 
of thousands of nearby galaxies in controlled surveys and probed their lumi- 
nous precursors out to redshift z ~6 - corresponding to a period only 1 Gyr 
after the Big Bang. Most importantly, a standard model has emerged which, 
through detailed numerical simulations, is capable of detailed predictions and 
interpretation of obscrvablcs. Many puzzles remain, as we will see, but the 
progress is truly impressive. 

This gives us confidence to begin addressing the final frontier in galaxy 
evolution: the earliest stellar systems and their influence on the intcrgalactic 
medium. When did the flrst substantial stellar systems begin to shine? Were 
they responsible for reionizing hydrogen in intcrgalactic space and what phys- 
ical processes occurring during these early times influenced the subsequent 
evolution of normal galaxies? 
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Let's begin by considering a crude measure of our recent progress. Figure 
1 shows the rapid pace of discovery in terms of the relative fraction of the 
refereed astronomical Uterature in two North American journals pertaining 
to studies of galaxy evolution and cosmology. These are cast alongside some 
milestones in the history of optical facilities and the provision of widely-used 
datasets. The figure raises the interesting question of whether more publica- 
tions in a given field means most of the key questions are being answered. 
Certainly, we can conclude that more researchers are being drawn to work in 
the area. But some might argue that new students should move into other, 
less well-developed, fields. Indeed, the progress in cosmology, in particular, is 
so rapid that some have raised the specter that the subject may soon reaching 
some form of natural conclusion (c.f. Horgan 1998). 
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Fig. 1. Fraction of the refereed astronomical literature in two North American jour- 
nals related to galaxy evolution and the cosmological parameters. The survey implies 
more than a doubling in fractional share over the past 15 years. Some possibly- 
associated milestones in the provision of unique facilities and datasets are marked 
(courtesy: J. Brinchmann). 

I believe, however, that the rapid growth in the share of publications is 
largely a reflection of new-found observational capabilities. We are witnessing 
an expansion of exploration which will most likely be followed with a more 
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detailed physical phase where we will be concerned with understanding how 
galaxies form and evolve. 

1.2 Observations Lead to Surprises 

It's worth emphasizing that many of the key features which define our current 
view of the Universe were either not anticipated by theory or initially rejected 
as unreasonable. Here is my personal short list of surprising observations which 
have shaped our view of the cosmos: 

1. The cosmic expansion discovered by Slipher and Hubble during the pe- 
riod 1917-1925 was not anticipated and took many years to be accepted. 
Despite the observational evidence and the prediction from General Rela- 
tivity for evolution in world models with gravity, Einstein maintained his 
preference for a static Universe until the early 1930's. 

2. The hot Big Bang picture received widespread support only in 1965 upon 
the discovery of the cosmic microwave background (Penzias & Wilson 
1965). Although many supported the hypothesis of a primeval atom, Hoyle 
and others considered an unchanging 'Steady State' universe to be a more 
natural solutuion. 

3. Dark matter was inferred from the motions of galaxies in clusters over 
seventy years ago (Zwicky 1933) but no satisfactory explanation of this 
puzzling problem was ever presented. The ubiquity of dark matter on 
galactic scales was realized much later (Rubin et al 1976). The dominant 
role that dark matter plays in structure formation only followed the recent 
observational evidence (Blumenthal et al 1984)0 

4. The cosmic acceleration discovered independently by two distant super- 
novae teams (Riess et al 1998, Perlmutter et al 1999) was a complete 
surprise (including to the observers, who set out to measure the deceler- 
ation). Although the cosmological constant, A, had been invoked many 
times in the past, the presence of dark energy was completely unforeseen. 

Given the observational opportunities continue to advance, it seems rea- 
sonable to suppose further surprises may follow! 

1.3 Recent Observational Milestones 

Next, it's helpful to examine a few of the most significant observational 
achievements in cosmology and structure formation over the past ~15 years. 
Each provides the basis of knowledge from which we can move forward, elim- 
inating a range of uncertainty across a wide field of research. 



For an amusing musical history of the role of dark mat- 
ter in cosmology suitable for students of any age check out 
: //www-astronomy .mps . ohio-state . edu-^dhw/Silliness/silliness .htmll 
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The Rate of Local Expansion: Bubble's Constant 

The Hubble Space Telescope (HST) was partly launched to resolve the puz- 
zling dispute between various observers as regards to the value of Hubble's 
constant Hq, normally quoted in kms sec~^ Mpc~^, or as h, the value in units 
of 100 kms sec""'^ Mpc"^. During the planning phases, a number of scientific 
key projects were defined and proposals invited for their execution. 

A very thorough account of the impasse reached by earlier ground-based 
observers in the 1970's and early 1980's can be found in Rowan-Robinson 
(1985) who reviewed the field and concluded a compromise of 67 ± 15 kms 
sec~^ Mpc~^, surprisingly close to the presently-accepted value. Figure 2 
nicely illustrates the confused situation. 
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Fig. 2. Various values of Hubble's constant in units of kms sec~^ Mpc~^ plotted as 
a function of the date of publication. Labels refer to estimates by Sandage & Tam- 
mann, de Vaucouleurs, van den Bergh and their respective collaborators. Estimates 
from the HST Key Project group (Freedman et al 2001) are labeled KP. From an ini- 
tial range spanning 50< Ho <100, a gradual convergence to the presently-accepted 
value is apparent. (Plot compiled and kindly made available by J. Huchra) 



Figure 3 shows the two stage 'step-ladder' technique used by Freedman et 
al (2001) who claim a final value of 67 ± 15 kms sec~^ Mpc~^. 'Primary' dis- 
tances were estimated to a set of nearby galaxies via the measured brightness 
and periods of luminous Cepheid variable stars located using HST's WFPC-2 
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imager. Over the distance range across which such individual stars can be 
seen (<25 Mpc), the leverage on Hq is limited and seriously affected by the 
peculiar motions of the individual galaxies. At ~20 Mpc, the smooth cosmic 
expansion would give Vexp —1400 kms sec~^ and a 10% error in Hq would 
provide a comparable contribution, at this distance, to the typical peculiar 
motions of galaxies of Vpec —50-100 kms sec~^. Accordingly, a secondary dis- 
tance scale was established for spirals to 400 Mpc distance using the empirical 
relationship first demonstrated by TuUy & Fisher (1977) between the /-band 
luminosity and rotational velocity. At 400 Mpc, the effect of Vpec is negligible 
and the leverage on Hq is excellent. Independent distance estimators utiliz- 
ing supernovae and elliptical galaxies were used to verify possible systematic 
errors. 




Distance (MpG) DlsUnce (Mpc) 



Fig. 3. Two step approach to measuring Hubble's constant Ho - the local expansion 
rate (Freedman et al 2001). (Left) Distances to nearby galaxies within 25 Mpc 
were obtained by locating and monitoring Cepheid variables using HST's WFPC-2 
camera; the leverage on Hq is modest over such small distances and affected seriously 
by peculiar motions. (Right) Extension of the distance-velocity relation to 400 Mpc 
using the I-band TuUy-Fisher relation and other techniques. The absolute scale has 
been calibrated using the local Cepheid scale. 



Cosmic Microwave Background: Thermal Origin and Spatial 
Flatness 

The second significant milestone of the last 15 years is the improved under- 
standing of the cosmic microwave background (CMB) radiation, commencing 
with the precise black body nature of its spectrum (Mather et al 1990) in- 
dicative of its thermal origin as a remnant of the cosmic fireball, and the sub- 
sequent detection of fluctuations (Smoot et al 1992), both realized with the 
COBE satellite data. The improved angular resolution of later ground-based 
and balloon-borne experiments led to the isolation of the acoustic horizon 
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scale at the epoch of recombination (de Bernadis et al 2000, Hanany et al 
2000). Subsequent improved measures of the angular power spectrum by the 
Wilkinson Microwave Anisotropy Probe (WMAP, Spergel et al 2003, 2006) 
have refined these early observations. The location of the primary peak in the 
angular power spectrum at a multiple moment I ~200 (corresponding to a 
physical angular scale of ~1 degree) provides an important constraint on the 
total energy density Qtot and hence spatial curvature. 

The derivation of spatial curvature from the angular location of the first 
acoustic (or 'Doppler') peak. Or, is not completely independent of other cos- 
mological parameters. There are dependences on the scale factor via Hq and 
the contribution of gravitating matter 1?m, viz: 

where h is Hq in units of 100 kms sec~^ Mpc~^. 

However, in the latest WMAP analysis, combining with distant supernovae 
data, space is flat to within 1%. 

Clustering of Galaxies: Gravitational Instability 

Galaxies represent the most direct tracer of the rich tapestry of structure in 
the local Universe. The 1970's saw a concerted effort to introduce a formalism 
for describing and interpreting their statistical distribution through angular 
and spatial two point correlation functions (Peebles 1980). This, in turn, led 
to an observational revolution in cataloging their distribution, first in 2-D 
from panoramic photographic surveys aided by precise measuring machines, 
and later in 3-D from multi-object spectroscopic redshift surveys. 

The angular correlation function w{8) represents the excess probability 5 P 
of finding a pair of galaxies separated by an angular separation (degrees). 

In a catalog averaging N galaxies per square degree, the probability of 
finding a pair separated by can be written: 

6 P = N[l + w{e)\S Q 

where 5 O is the solid angle of the counting bin, (i.e. 9 to 9 + 69). 
The corresponding spatial equivalent, ^(r) in a catalog of mean density p 
per Mpc^ is thus: 

SP^ p[l+^{r)]SV 

One can be statistically linked to the other if the overall redshift distribu- 
tion of the sources is available. 

Figure 4 shows a pioneering detection of the angular correlation function 
w{0) for the Cambridge APM Galaxy Catalog (Maddox et al 1990). This was 
one of the first well-constructed panoramic 2-D catalogs from which the large 
scale; nature of the galaxy distribution could be discerned. A power law form 
is evident: 
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where, for example, 9 is measured in degrees. The amplitude A decreases 
with increasing depth due to both increased projection from physically- 
uncorrelated pairs and the smaller projected physical scale for a given angle. 




Fig. 4. Angular correlation function for the APM galaxy catalog - a photographic 
survey of the southern sky (Maddox et aJ 1990) - partitioned according to limiting 
magnitude (left). The amplitude of the clustering decreases with increasing depth 
due to an increase in the number of uncorrelated pairs and a smaller projected 
physical scale for a given angle. These effects can be corrected in order to produce 
a high signal/noise function scaled to a fixed depth clearly illustrating a universal 
power law form over nearly 3 dex (right). 



Highly-multiplexed spectrographs such as the 2 degree field instrument on 
the Anglo-Australian Telescope (Colless et al 2001) and the Sloan Digital Sky 
Survey (York et al 2001) have led to the equivalent progress in 3-D surveys 
(Figure 5). In the early precursors to these grand surveys, the 3-D equivalent 
of the angular correlation function, was also found to be a power law: 

ro 

where (Mpc) is a valuable clustering scale length for the population. 

As the surveys became more substantial, the power spectrum P{k) has 
become the preferred analysis tool because its form can be readily predicted 
for various dark matter models. For a given density field p(x), the fluctuation 



Observations of the High Redshift Universe 



9 



over the mean is d = p / p and for a given wavenumber k, the power spectrum 
becomes: 

P{k) =< \Sl\ >= J £_{r)exp{i^^.r)d^r 

The final power spectrum for the completed 2dF survey is shown in Figure 
6 (Cole et al 2005) and is in remarkably good agreement with that predicted 
for a cold dark matter spectrum consistent with that which reproduces the 
CMB angular fluctuations. 




Fig. 5. Galaxy distribution from the completed 2dF redshift survey (CoUess et al 
2001). 



Dark Matter and Gravitational Instability 

We have already mentioned the ubiquity of dark matter on both cluster and 
galactic scales. The former was recognized as early as the 1930's from the high 
line of sight velocity dispersion aios of galaxies in the Coma cluster (Zwicky 
1933). Assuming simple virial equilibrium and isotropically-arranged galaxy 
orbits, the cluster mass contained with some physical scale Rd is: 

M = 3 < al, > Rei/G 

which far exceeds that estimated from the stellar populations in the cluster 
galaxies. High cluster masses can also be confirmed completely independently 
from gravitational lensing where a background source is distorted to produce 
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Fig. 6. Power spectrum from the completed 2dF redshift survey (Cole et al 2005). 
Solid lines refer to the input power spectrum for a dark matter model with the tab- 
ulated parameters and that convolved with the geometric 'window function' which 
affects the observed shape on large scales. 



a 'giant arc' - in effect a partial or incomplete 'Einstein ring' whose diameter 
0E for a concentrated mass M approximates: 



and D = Di , / Dds where the subscripts s and I refer to angular diam- 
eters distances of the background source and lens respectively. 

On galactic scales, extended rotation curves of gaseous emission lines in 
spirals (see review by Rubin 2000) can trace the mass distribution on the 
assumption of circular orbits, viz: 

GM{< R) _ 

Flat rotation curves {V ^constant) thus imply A/(< R) oc R. Together 
with arguments based on the question on the stability of flattened disks (Os- 
triker & Peebles 1973), such observations were critical to the notion that all 
spiral galaxies are embedded in dark extensive 'halos'. 

The evidence for halos around local elliptical galaxies is less convincing 
largely because there are no suitable tracers of the gravitational potential on 
the necessary scales (see Gerhard et al 2001). However, by combining grav- 
itational lensing with stellar dynamics for intermediate redshift ellipticals, 
Koopmans & Treu (2003) and Treu et al (2006) have mapped the projected 
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dark matter distribution and show it to be closely fit by an isothermal profile 
p{r) oc r~^. 

The presence of dark matter can also be deduced statistically from the dis- 
tortion of the galaxy distribution viewed in redshift space, for example in the 
2dF survey (Peacock ct al 2001). The original idea was discussed by Kaiser et 
al (1987). The spatial correlation function ^(r) is split into its two orthogonal 
components, ^(cr, tt) where a represents the projected separation perpendicu- 
lar to the line of sight (unaffected by peculiar motions) and tt is the separation 
along the line of sight (inferred from the velocities and hence used to mea- 
sure the effect). The distortion of ^(ct, tt) in the tt direction can be measured 
on various scales and used to estimate the line of sight velocity dispersion 
of pairs of galaxies and hence their mutual gravitational field. Depending on 
the extent to which galaxies are biased tracers of the density field, such tests 
indicate i7M=0.25. 

On the largest scales, weak gravitational lensing can trace the overall distri- 
bution and dark matter content of the Universe (Blandford & Narayan 1992, 
Refrcgicr 2003). Recent surveys are consistent with these estimates (Hoekstra 
et al 2005). 



Dark Energy and Cosmic Acceleration 

Prior to the 1980's observational cosmologists were obsessed with two empiri- 
cal quantities though to govern the cosmic expansion history - R{t): Bubble's 
constant Hq = dR/dt and a second derivative, the deceleration parameter qq, 
which would indicate the fate of the expansion: 

_ d^R/dt^ 
^° ~ (dR/dt)^ 

In the presence only of gravitating matter, Priedmann cosmologies indicate 
/2m = 2(7o- The distant supernovac searches were begun in the expectation of 
measuring qo independently of f2M and verifying a low density Universe. 

As we have discussed, Type la supernovae (SNe) were found to be fainter 
at a given recessional velocity than expected in a Universe with a low mass 
density; Figure 7 illustrates the effect for the latest results from the Canada- 
France SN Legacy Survey (Astier et al 2006). In fact the results cannot be 
explained even in a Universe with no gravitating matter! A formal fit for qq 
indicates a negative value corresponding to a cosmic acceleration. 

Acceleration is permitted in Priedmann models with a non-zero cosmolog- 
ical constant A. In general (Caroll et al 1992): 

where fi^ = yl/87r G is the energy density associated with the cosmological 
constant. 
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Fig. 7. Hubble diagram (distance-redshift relation) for calibrated Type la super- 
novae from the first year data taken by the Canada France Supernova Legacy Survey 
(Astier et al 2006). Curves indicate the relation expected for a high density Universe 
without a cosmological constant and that for the concordance cosmology (see text) 
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The appeal of resurrecting the cosmological constant is not only its ability 
to explain the supernova data but also the spatial flatness in the acoustic peak 
in the CMB through the combined energy densities /2m + f^A - the so-called 
Concordance Model (Ostriker & Steinhardt 1995, Bahcall et al 1999). 

Hovifcvcr, the observed acceleration raises many puzzles. The absolute value 
of the cosmological constant cannot be understood in terms of physical de- 
scriptions of the vacuum energy density, and the fact that i?M — implies 
the accelerating phase began relatively recently (at a redshift of Z2±0.7). Al- 
ternative physical descriptions of the phenomenon (termed 'dark energy') are 
thus being sought which can be generalized by imagining the vacuum obeys an 
equation of state where the negative pressure p relates to the energy density 
p via an index w, 

p = w p 

in which case the dependence on the scale factor R goes as 

p OC J?-3(l+-) 

The case •w=-l would thus correspond to a constant term equivalent to 
the cosmological constant, but in principle any w; <-l/3 would produce an 
acceleration and conceivably w is itself a function of time. The current SNLS 
data indicate w=-1.023 ± 0.09 and combining with the WMAP data does not 
significantly improve this constraint. 

1.4 Concordance Cosmology: Why is such a curious model 
acceptable? 

According to the latest WMAP results (Spergel et al 2006) and the analysis 
which draws upon the progress reviewed above (the HST Hubble constant Key 
Project, the large 2dF and SDSS redshift surveys, the CFHT supernova survey 
and the first weak gravitational lensing constraints), we live in a universe with 
the constituents listed in Table 1. 



Table 1. Cosmic Constituents 



Total Matter 




0.24 


± 


0.03 


Baryonic Matter 


Ob 


0.042 


± 


0.004 


Dark Energy 


Oa 


0.73 


± 


0.04 



Given only one of the 3 ingredient is physically understood it may be 

reasonably questioned why cosmologists are triumphant about having reached 
the era of 'precision cosmology' ! Surely we should not confuse measurement 
with understanding? 
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The underlying reasons are two-fold. Firstly, many independent probes 
(rcdshift surveys. CMB fluctuations and Icnsing) indicate the low matter den- 
sity. Two independent probes not discussed (primordial nucleosynthesis and 
CMB fluctuations) support the baryon fraction. Finally, given spatial flatness, 
even if the supernovae data were discarded, we would deduce the non-zero dark 
energy from the above results alone. 

Secondly, the above parameters reconcile the growth of structure from the 
CMB to the local redshift surveys in exquisite detail. Numerical simulations 
based on 10^° particles (e.g. Springel et al 2005) have reached the stage where 
they can predict the non-linear growth of the dark matter distribution at 
various epochs over a dynamic range of 3-4 dcx in physical scales. Although 
some input physics is needed to predict the local galaxy distribution, the 
agreement for the concordance model (often termed ylCDM) is impressive. In 
short, a low mass density and non-zero A both seem necessary to explain the 
present abundance and mass distribution of galaxies. Any deviation would 
either lead to too much or too little structure. 

This docs not mean that the scorecard for /ICDM should be considered 
perfect at this stage. As discussed, we have little idea what the dark matter or 
dark energy might be. Moreover, there are numerous difiiculties in reconcil- 
ing the distribution of dark matter with observations on galactic and cluster 
scales and frequent challenges that the mass assembly history of galaxies is 
inconsistent with the slow hierarchical growth expected in a yl-dominated Uni- 
verse. However, as wc will see in later lectures, most of these problems relate 
to applications in environments where dark matter co-exists with baryons. 
Understanding how to incorporate baryons into the very detailed simulations 
now possible is an active area where interplay with observations is essential. 
It is helpful to view this interplay as a partnership between theory and obser- 
vation rather than the oft-quoted 'battle' whereby observers challenge or call 
into question the basic principles. 

1.5 Lecture Summary 

I have spent my first lecture discussing largely cosmological progress and the 
impressive role that observations have played in delivering rapid progress. 

All the useful cosmological functions - e.g. time, distance and comoving 
volume versus redshift, are now known to high accuracy which is tremendously 
beneficial for our task in understanding the first galaxies and stars. I emphasize 
this because oven a decade ago, none of the physical constants were known 
well enough for us to be sure, for example, the cosmic age corresponding to a 
particular redshift. 

I have justified A CDM as an acceptable standard model, despite the mi- 
known nature of its two dominant constituents, partly because there is a con- 
cordance in the parameters when viewed from various observational probes, 
and partly because of the impressive agreement with the distribution of galax- 
ies on various scales in the present Universe. 
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Connecting the dark matter distribution to the observed properties of 
galaxies requires additional physics relating to how baryons cool and form 
stars in dark matter halos. Detailed observations are necessary to 'tune' the 
models so these additional components can be understood. 

All of this will be crucial if we are correctly predict and interpret signals 
from the first objects. 
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2 Galaxies &c the Hubble Sequence 

2.1 Introduction: Changing Paradigms of GalcLxy Formation 

Wc now turn to the interesting history of how our views of galaxy formation 
have changed over the past 20-30 years. It is convenient to break this into 3 
eras 

1. The classical era (pre-1985) as articulated for example in the influential 

articles by Beatrice Tinsley and others. Galaxies were thought to evolve in 
isolation with their present-day properties governed largely by one func- 
tion - the time-dependent star formation rate tp{t). Ellipticals suffered a 
prompt conversion of gas into stars, whereas spirals were permitted a more 
gradual consumption rate leading to a near-constant star formation rate 
with time. 

2. The dark matter-based era (1985-): in hierarchical models of structure 
formation involving gravitational instability, the ubiquity of dark matter 
halos means that merger driven assembly is a key feature. If mergers 
redistribute angular momentum, galaxy morphologies arc transformed. 

3. Understanding feedback and the environment (1995-): In the most re- 
cent work, the evolution of the morphology-density relation (Dressier et 
al 1997) and the dependence of the assembly history on galactic mass 
('downsizing', Cowie et al 1996) have emphasized that star formation is 
regulated by processes other than gas cooling and infall associated with 
DM- driven mergers. 

2.2 Galaxy Morphology - Valuable Tool or Not? 

In the early years, astronomers placed great stock on understanding the ori- 
gin of the morphological distribution of galaxies, sometimes referred to as the 
Hubble sequence (Hubble 1936). Despite this simple categorization 70 years 
ago, the scheme is evidently still in common use. In its support, Sandage (e.g. 
2005) has commented on this classification scheme as describing 'a true order 
among the galaxies, not one imposed by the classifier'. However, many con- 
temporary modelers and observers have paid scant attention to morphology 
and placed more emphasis on understanding stellar population differences. 
What value should we place on accurately measuring and reproducing the 
morphological distribution? 

The utility of Hubble's scheme, at least for local galaxies, lies in its ability 
to distinguish dynamically distinct structures - spirals and SOs arc rotating 
stellar disks, whereas luminous spheroids are pressure-supported ellipsoidal or 
triaxial systems with anisotropic velocity fields. This contains key information 
on the degree of dissipation in their formation (Fall & Efstathiou 1980). 

There are also physical variables that seem to underpin the sequence, 
including (i) gas content and color which relate to the ratio of the current to 
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past average star formation ratio ip{ta)/ ip (Figure 8) and (ii) inner structures 
including the bulge-to-disk ratio. Various modelers (Baugh et al 1996) have 
argued that the bulge-to-disk ratio is closely linked to the merger history and 
attempted to reproduce the present distribution as a key test of hierarchical 
assembly. 
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Fig. 8. A succinct summary of the classical view of galaxy formation (pre-1985): 
(Left) The monotonia distribution of Hubble sequence galax;ies in the U—VvsV—K 
color plane (Aaaronson 1978). (Right) A simple model which reproduces this trend 
by changing only the ratio of the current to past average star formation rate (Struck- 
Marcell & Tinsley 1978). Galaxies with constant star formation permanently occupy 
the top left (blue) corner; galaxies with an initial burst rapidly evolve to the bottom 
right (red) corner. 



Much effort has been invested in attempting to classify galaxies at high 
redshift, both visually and with automated algorithms. This is a challenging 
task because the precise appearance of diagnostic features such as spiral arms 
and the bulge/disk ratio depends on the rest- wavelength of the observations. 

An effect termed the 'niorphological k-correction' can thus shift galaxies to 
apparently later types as the redshift increases for observations conducted 
in a fixed band. A further limitation, which works in the opposite sense, is 
surface brightness dimming, which proceeds as oc (1 + z)*, rendering disks 
less prominent at high redshift and shifting some galaxies to apparent earlier 
types. 

The most significant achievements from this effort has been the realization 
that, despite the above quantitative uncertainties, faint star-forming galaxies 
are generally more irregular in their appearance than in local samples (Glaze- 
brook et al 1995, Driver et al 1995). Moreover, HST images suggest on-going 
mergers with an increasing frequency at high redshift (LeFevre et al 2000) al- 
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though quantitative estimates of the merging fraction as a function of redshift 
remain uncertain (sec Bundy ct al 2004). 

The idea that morphology is driven by mergers took some time for the ob- 
servational community to accept. Numerical simulations by Toomre & Toomre 
(1972) provided the initial theoretical inspiration, but the observational ev- 
idence supporting the notion that spheroidal galaxies were simple collapsed 
systems containing old stars was strong (Bower et al 1992). Tell-tale signs 
of mergers in local ellipticals include the discovery of orbital shells (Malin 
& Carter 1980) and multiple cores revealed only with 2-D dynamical studies 
(Davies et al 2001). 

2.3 Semi-Analytical Modeling 

As discussed by the other course lecturers and briefly in §1, our ability to follow 
the distribution of dark matter and its growth in numerical simulations is well- 
advanced (e.g. Springel et al 2005). The same cannot be said of understanding 
how the baryons destined, in part, to become stars are allocated to each DM 
halo. This remains the key issue in interfacing theory to observations. 

Progress has occurred in two stages - according to the eras discussed 
in §2.1. Semi analytic codes were first developed in the 1990's to introduce 
baryons into DM n-body simulations using prescriptive methods for star for- 
mation, feedback and morphological assembly (Figure 9). These codes were 
initially motivated to demonstrate that the emerging DM paradigm was con- 
sistent with the abundance of observational data (Kauffmann et al 1993, 
Somerville & Primack 1999, Cole et al 2000). Prior to development of these 
codes, evolutionary predictions were based almost entirely on the 'classical' 
viewpoint with stellar population modeling based on variations in the star 
formation history 'tp{t) for galaxies evolving in isolation e.g. Bruzual (1980). 

Initially these feedback prescriptions were adjusted to match observables 
such as the luminosity function (whose specific details we will address below), 
as well as specific attributes of various surveys (counts, redshift distributions, 
colors and morphologies). In the recent versions, more elaborate physically- 
based models for feedback processes are being considered (e.g. Croton et al 
2006) 

The observational community was fairly skeptical of the predictions from 
the first semi- analytical models since it was argued that the parameter space 
implied by Figure 9 enabled considerable freedom even for a fixed primor- 
dial fluctuation spectrum and cosmological model. Moreover, where different 
codes could be compared, considerably different predictions emerged (Benson 
et al 2002). Only as the observational data has moved from colors and star 
formation rates to physical variables more closely related to galaxy assem- 
bly (such as stellar masses) have the limitations of the early semi-analytical 
models been exposed. 
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Fig. 9. Schematic of the ingredients inserted into a semi-analytical model (from 
Cole et al 2000). Solid lines refer to mass transfer, dashed lines to the transfer of 
metals according to different compositions Z. Gas cooling (M) and star formation 
(?/;) is inhibited by the effect of supernovae (/3). Stars return some fraction of their 
mass to the interstellar medium (i?) and to the hot gas phase (e) according to a 
metal yield p. 

2.4 A Test Case: The Galaxy Luminosity Function 

One of the most straightforward and fundamental predictions a theory of 
galaxy formation can make is the present distribution of galaxy luminosities 
- the luminosity function (LF) ^(L) whose units are normally per comoving 
Mpc^B 

As the contribution of a given luminosity bin dL to the integrated luminos- 
ity density per unit volume is cx ^(V) L dL, an elementary calculation shows 
that all luminosity functions (be they for stars, galaxies or QSOs) must have a 
bend at some characteristic luminosity, otherwise they would yield an infinite 
total luminosity (see Felton 1977 for a cogent early discussion of the signifi- 
cance and intricacies of the LF). Recognizing this, Schechter (1976) proposed 
the product of a power law and an exponential as an appropriate analytic 
representation of the LF, viz: 

<|>(i)^^^*(iL)-"exp(-i^)^ 

where ^* is the overall normalization corresponding to the volume density 
at the turn-over (or characteristic) luminosity L* , and a is the faint end slope 
which governs the relative abundance of faint and luminous galaxies. 

The total abundance of galaxies per unit volume is then: 



^ If Hubble's constant is not assumed, it is quoted in units of ft ^ Mpc 
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Ntot = J ^(i) dL = <P* r{a + 1) 
and the total luminosity density is: 

PL = J ${L)LdL = r{a + 2) 

where F is the incomplete gamma function which can be found tabulated 
in most books with integral tables (e.g. Gradshteyn & Ryzhik 2000). Note 
that Ntot diverges if a <-l, whereas pL diverges only if a <-2. 

Recent comprehensive surveys by the 2dF team (Norberg et al 2002) and 
by the Sloan Digital Sky Survey (Blanton et al 2001) have provided definitive 
values for the LF in various bands. Encouragingly, when allowance is made 
for the various photometric techniques, the two surveys are in excellent agree- 
ment. Figure 10 shows the Schcchter function is a reasonably good (but not 
perfect) fit to the 2dF data limited at apparent magnitude hj <19.7. Moreover 
there is no significant difference between the LF derived independently for the 
two Galactic hemispheres. The slight excess of intrinsically faint galaxies in 
the northern cap is attributable to a local inhomogeneity in the nearby Virgo 
supcrcluster. 

Fundamental though this function is, despite ten years of semi-analytical 
modeling, reproducing its form has proved a formidable challenge (as discussed 
by Benson et al 2003, Croton et al 2006 and de Lucia et al 2006). Early 
predictions also failed to reproduce the color distribution along the LF. The 
halo mass distribution does not share the sharp bend at L* and too much star 
formation activity is retained in massive galaxies. These early predictions 
produced too many luminous blue galaxies and too many faint red galaxies 
(Bower et al 2006). 

As a result, more specific feedback recipes have been created to resolve 
this discrepancy. Several physical processes have been invoked to regulate 
star formation as a function of mass, viz: 

• Reionization feedback: radiative heating from the first stellar systems at 
high redshift which increases the Jeans mass, inhibiting the early formation 
of low mass systems, 

• Supernova feedback: this was considered in the early semi-analytical mod- 
els but is now more precisely implemented so as to re-heat the interstellar 
medium, heat the halo gas or even eject the gas altogether from low mass 
systems, 

• Feedback from active galactic nuclei: the least well-understood process 
with various modes postulated to transfer energy from an active nucleus 
to the halo gas. 

Benson et al (2003) and Croton et al (2006) illustrate the effects of these 
more detailed prescriptions for these feedback modes on the predicted LF and 
find that supernova and reionization feedback largely reduce the excess of 
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Fig. 10. Rcst-framc bj luminosity function from the 2dF galaxy redshift survey 
(Norberg et al 2002). (Top) a comparison of results across the northern and southern 
Galactic caps; there is only a marginal difference in the abundance of intrinsically 
faint galaxies. (Bottom) combined from both hemispheres indicating the Schechter 
function is a remarkably good fit except at the extreme ends of the luminosity 
distribution. 

intrinsically faint galaxies but, on grounds of energetics, only AGN can inhibit 

star formation and continued growth in massive galaxies. There remains an 
excess at the very faint end (Fig. 11). 

2.5 The Role of the Environment 

In addition to recognizing that more elaborate modes of feedback need to be 
incorporated in theoretical models, the key role of the environment has also 
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Fig. 11. The effect of various forms of feedback (dots) on the resultant shape of the 
blue 2dF galaxy luminosity function (Norberg et al 2002). From left to right: no feed- 
back, supernova feedback only, supernova, reionization & AGN feedback (Courtesy: 
Darren Croton). 



emerged as an additional feature which can truncate star formation and alter 
galaxy morphologies. 

The preponderance of elliptical and SO galaxies in rich clusters was no- 
ticed in the 1930's but the first quantitative study of this effect was that of 
Dressier (1980) who correlated the fraction of galaxies of a given morphology 
T above some fixed luminosity with the projected galaxy density, S, measured 
in galaxies Mpc~^. 

The local T ^ S relation was used to justify two rather different possi- 
bilities. In the first, the nature hypothesis, those galaxies which formed in 
high density peaks at early times were presumed to have consumed their gas 
efficiently, perhaps in a single burst of star formation. Galaxies in lower den- 
sity environments continued to accrete gas and thus show later star formation 
and disk-like morphologies. In short, segregation was established at birth and 
the present relation simply represents different ways in which galaxies formed 
according to the density of the environment at the time of formation. In the 
second, the nurture hypothesis, galaxies are transformed at later times from 
spirals into spheroidals by environmentally-induced processes. 

Work in the late 1990's, using morphologies determined using Hubble 
Space Telescope, confirmed a surprisingly rapid evolution in the T—E relation 
over 0< z <0.5 (Couch et al 1998, Dressier et al 1997) strongly supporting 
environmentally-driven evolution along the lines of the nurture hypothesis. 
Impressive Hubble images of dense clusters at quite modest redshifts [z ~0.3- 
0.4) showed an abundance of spirals in their cores whereas few or none exist 
in similar environs today. 

What physical processes drive this relation and how has the T — S relation 
evolved in quantitative detail? Recent work (Smith et al 2005, Postman et al 
2005, Figure 12) has revealed that the basic relation was in place at z ~1, 
but that the fraction fs+so of Es and SOs has doubled in dense environments 
since that time. Smith et al suggest that a continuous, density-dependent. 
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transformation of spirals into SOs would explain the overall trend. Treu et al 
(2003) likewise see a strong dependence of the fraction as a function of S (and 
to a lesser extent with cluster-centric radius) within a well-studied cluster at 
z ~0.4; they review the various physical mechanisms that may produce such 
a transformation. 
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Fig. 12. The fraction of observed E and SO galaxies down to a fixed rest-frame 
luminosity as a function of lookback time and projected density S from the study of 
Smith et al (2005). Although the morphology density relation was already in place 
at z ~1, there has been a continuous growth in the fraction subsequently, possibly 
as a result of the density-dependent transformation of spirals into SOs. 



Figure 12 encapsulates much of what we now know about the role of the 
environment on galaxy formation. The early development of the T — re- 
lation implies dense peaks in the dark matter distribution led to accelerated 
evolution in gas consumption and stellar evolution and this is not dissimilar to 
the nature hypothesis. However, the subsequent development of this relation 
since z ~1 reveals the importance of environmentally-driven morphological 
transformations. 



2.6 The Importance of High Redshift Data 

This glimpse of evolving galaxy populations to z ~1 has emphasized the 
important role of high redshift data. In the case of the local morphology- 
density relation, Hubble morphologies of galaxies in distant clusters have given 
us a clear view of an evolving relationship, partly driven by environmental 
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processes. Indeed, the data seems to confirm the nurture hypothesis for the 
origin of the morphology-density relation. 

Although we can place important constraints on the past star formation 
history from detailed studies of nearby galaxies, as the standard model now 
needs several additional ingredients (e.g. feedback) to reproduce even the most 
basic local properties such as the luminosity function (Figure 11), data at 
significant look-back times becomes an essential way to test the validity of 
these more elaborate models. 

Starting in the mid-1990's, largely by virtue of the arrival of the Keck 
telescopes - the first of the new generation of 8-lOm class optical/infrared 
telescopes - and the refurbishment of the Hubble Space Telescope, there has 
been an explosion of new data on high redshift galaxies. 

It is helpful at this stage to introduce three broad classes of distant ob- 
jects which will feature significantly in the next few lectures. Each gives a 
complementary view of the galaxy population at high redshift and illustrates 
the challenge of developing a unified vision of galaxy evolution. 
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Fig. 13. Location and verification of the Lyman break population in the redshift 
range 2.7< z <3.4 (Steidel et al 2003). (Left) UGR color-color plane for a single field; 
green and yellow shading refers to variants on the color selection of high redshift 
candidates. Contaminating Galactic stars cut across the lower right corner of this 
selection; those confirmed spectroscopically are marked in red. (Right) Coadded 
rest-frame Keck spectra for samples of typically 200 Lyman break galaxies binned 
according to the strength of Lyman a emission. 
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• Lyman-break galaxies: color-selected luminous star forming galaxies at 
z >2. First located spectroscopically by Steidel et al (1996, 1999a, 1999b, 
2003), these sources are selected by virtue of the increased opacity short- 
ward of the Lyman limit (A =912A ) arising from the combined effect 
of neutral hydrogen in hot stellar atmospheres, the interstellar gas and 
the intergalactic medium. When redshifted beyond z~2, the characteristic 
'drop out' in the Lyman continuum moves into the optical (Figure 13). We 
will review the detailed properties of this, the most well-studied, distant 
galaxy population over 2< z <5 in subsequent lectures. 
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Fig. 14. Redshift distribution of 73 radio-identified SCUBA sources from Chapman 
et al (2005). To illustrate the possible bias arising from the necessary condition of 
a radio position for a Keck redshift, the solid curve represents a model prediction 
for the entire >5mJy sub- mm population. The sub-mm population does not seem 
to extend significantly beyond z~4 and has a median redshift of z=2.2. 

• Sub-millimeter star forming sources: The SCUBA 850/xm array on the 15 
meter James Clerk Maxwell Telescope and other sub-mm imaging devices 
have also been used to locate distant star forming galaxies (Small et al 
1997, Hughes et al 1998). In this case, emission is detected from dust, 
heated either by vigorous star formation or an active nucleus. Remarkably, 
their visibility does not fall off significantly with redshift because they are 
detected in the Raylcigh- Jeans tail of the dust blackbody spectrum (Blain 
et al 2002). 
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Progress in understanding the role and nature of this population has been 
slower because sub-mm sources are often not visible at optical and near- 
infrared wavelengths (due to obscuration) and the positional accuracy of 
the sub-mm arrays is too coarse for follow-up spectroscopy. The impor- 
tance of sub-mm sources lies in the fact that they contribute significantly 
to the star formation rate at high redshift. Regardless of their redshift, the 
source density at faint limits is 1000 times higher than a no-evolution pre- 
diction based on the local abundance of dusty IRAS sources. For several 
years the key issue was to nail the redshift distribution. 
Progress has been made by securing accurate positions using radio interfer- 
ometers such as the VLA (Fraycr ct al 2000). About 70% of those brighter 
than 5 mJy have VLA detections and spectroscopic redshift have now 
been determined for a significant fraction of this population (Chapman et 
al 2003, 2005, Figure 14). 




Wavelength (rim) 



Fig. 15. Observed spectral energy distribution of distant red sources selected with 
a red J — K color superimposed on model spectra (Pranx et al 2003) . Although some 

sources reveal modest star formation and can be spcctroscopically confirmed to lie 
above z ~2, others (such as the lower two examples) appear to be passively-evolving 
with no active star formation. 

• Passively-Evolving Sources: The Lyman-break and sub-mm sources are 
largely star-forming galaxies. The arrival of panoramic near-infrared cam- 
eras has opened the possibility of locating quiescent sources that are no 
longer forming stars. Such sources would not normally be detected via the 
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other techniques and so understanding their contribution to the integrated 

stellar mass at, say, z ~2, is very important. 

The nomenclature here is confusing with intrinsically red sources being 
termed 'extremely red objects (EROs)' or 'distant red galaxies (DRGs)' 
with no agreed selection criteria (sec McCarthy 2004 for a review). When 
star formation is complete, stellar evolution continues in a passive sense 
with main sequence dimming; the galaxy fades and becomes redder. 
Of particular interest arc the most distant examples which co-cxist along- 
side the sub-mm and Lyman-break galaxies, i.e. at z >2 selected according 
to their infrared J — K color (van Dokkum et al 2003, Figure 15). 

2.7 Lecture Summary 

We have seen in this brief tour that galaxy formation is a process involving 
gravitational instability driven by the hierarchical assembly of dark matter 
halos; this component we understand well. However, additional complexities 
arise from star formation, dynamical interactions and mergers, environmental 
processes and various forms of feedback which serve to regulate how star 
formation continues as galaxies grow in mass. 

Theorists have attempted to deal with this complexity by augmenting the 
highly-successful numerical (DM-only) simulations with semi-analytic tools 
for incorporating these complexities. As the datasets have improved so it is 
now possible to consider 'fine-tuning' these semi-analytical ingredients. Ah 
initio modeling is never likely to be practical. 

I think it fair to say that many observers have philosophical reservations 
about this 'fine-tuning' process in the sense that although it may be possible 
to reach closure on models and data, we seek a deeper understanding of the 
physical reality of many of the ingredients. This is particularly the case for 
feedback processes. Fortunately, high redshift data forces this reality check as 
it gives us a direct measure of the galaxy assembly history which will be the 
next topics we discuss. 

As a way of illustrating the importance of high redshift data, I have in- 
troduced three very different populations of galaxies each largely lying in 
the redshift range 2< z <4. When these were independently discovered, it 
was (quite reasonably) c;laimed by their discoverers that their category repre- 
sented a major, if not the most significant, component of the distant galaxy 
population. We now realize that UV-selected, sub-mm selected and non-star 
forming galaxies each provide a complementary view of the complex history 
of galaxy assembly and the challenge is to complete the 'jig-saw' from these 
populations. 



28 Richard S Ellis 



3 Cosmic Star Formation Histories 

3.1 When Did Galaxies Form? Secirches for Primeval GalcLxies 

The question of the appearance of an early forming galaxy goes back to the 
1960's. Partride & Peebles (1967) imagined the free-fall collapse of a 700 L* 
system at z ~10 and predicted a diffuse large object with possible Lyman a 
emission. Meier (1976) considered primeval galaxies might be compact and 
intense emitters such as quasars. 

In the late 1970's and 1980's when the (then) new generation of 4 meter 
telescopes arrived, astronomers sought to discover the distinct era when galax- 
ies formed. Stellar synthesis models (Tinsley 1980, Bruzual 1980) suggested 
present-day passive systems (E/SOs) could have formed via a high redshift lu- 
minous initial burst. Placed at z ~2-3, sources of the same stellar mass would 
be readily detectable at quite modest magnitudes, B ~22-23, and provide an 
excess population of blue galaxies. 

In reality, the (now well-studied) excess of faint blue galaxies over locally- 
based predictions is understood to be primarily a phenomenon associated with 
a gradual increase in star formation over < 2; <1 rather than one due to a 
distinct new population of intensely luminous sources at high redshift (Koo 
& Kron 1992, Elhs 1997). Moreover, dedicated searches for suitably intense 
Lyman a emitters were largely unsuccessful. Pritchet (1994) comprehensively 
reviews a decade of searching. 

Our thinking about primeval galaxies changed in two respects in the late 
1980's. Foremost, synthesis models such as those developed by Tinsley and 
Bruzual assumed isolated systems; dark matter-based models emphasized the 
gradual assembly of massive galaxies. This change meant that, at z ~2-3, 
the abundance of massive galaxies should be much reduced. Secondly, the 
flux limits searched for primeval galaxies were optimistically bright; we slowly 
realized the more formidable challenge of finding these enigmatic sources. 

3.2 Local Inventory of Stars 

An important constraint on the past star formation history is the present-day 
stellar density. The former must, when integrated, yield the latter. Fukugita 
et al (1998) and Fukugita & Peebles (2004) have considered this important 
problem based on local survey data provided by the SDSS (Kauffmann et al 
2003) and 2dF (Cole et al 2001) redshift samples. 

The derivation of the integrated density of stars involves many assump- 
tions and steps but is based primarily on the local infrared (JC-band) lumi- 
nosity function of galaxies. The rest-frame K luminosity of a galaxy is a much 
more reliable proxy for its stellar mass than that at a shorter (e.g. optical) 
wavelength because its value is largely irrespective of the past star formation 
history - a point illustrated by Kauffmann & Chariot (1998, Figure 16). An- 
other way to phrase this is to say that the infrared mass/light ratio {M/Lk) 
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is fairly independent of the star formation history, so that the stellar mass can 
be derived from the observed ii'-band luminosity by a multiplicative factor. 
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Fig. 16. The robustness of the /f-band luminosity of a galaxy as a proxy for stellar 

mass (Kauffmann & Chariot 1998). The upper panel shows the K-hand apparent 
magnitude of a galaxy defined to have a fixed stellar mass of 10^'' Mq when placed 
at various redshifts. The different curves represent extreme variations in the way 
the stellar population was created. Whereas the B — K color is strongly dependent 
upon the star formation history, the K-banA luminosity is largely independent of it. 



In practice the mass/light ratio depends on the assumed distribution of 
stellar masses in a stellar population. The zero age or initial mass function is 
usually assumed to be some form of power law which can only be determined 
reliable for Galactic stellar populations, although constraints are possible for 
extragalactic populations from colors and nebular line emission (see reviews 
by Scalo 1986, Kennicutt 1998, Chabrier 2003) 

In its most frequently-used form the IMF is quoted in mass fraction per 
logarithmic mass bin: viz: 

^(logm) = — oc 
dlogm 



or, occasionally. 



30 Richard S Ellis 



£_{m) = — = —n—rT^Wogm) oc m " 
dm m(m 10) 

where x — a ~ \. 

In his classic derivation of the IMF, Salpeter (1955) determined a pure 
power law with x=1.35. More recently adopted IMFs are compared in Figure 
17. They differ primarily in how to restrict the low mass contribution, but 
there is also some dispute on the high mass slope (although the Salpeter value 
is supported by various observations of galaxy colors and Ha distributions, 
Kennicutt 1998). 
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Fig. 17. A comparison of popular stellar initial mass functions (courtesy of Ivan 
Baldry). 



The IMF has a direct influence on the assumed M/Lk (as discussed by 
Baldry & Glazebrook 2003, Chabrier, 2003 and Fukugita & Peebles, 2004) in 
a manner which depends on the age, composition and past star formation his- 
tory. The adopted mass/light ratio is then a crucial ingredient for computing 
both stellar masses (LectureLecture 4) and galaxy colors. 

BaldrjH has undertaken a very useful comparative study of the impact of 
various IMF assumptions using the PEGASE 2.0 stellar synthesis code for 
a population 10 Gyr old with solar metallicity, integrating between stellar 



^ lLttp://www.astro.livjm.ac.uk~ikb/research/imf-use-in-cosmology.litml 
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Table 2. K-band Stellar Mass/Light Ratios 



Source 


Stars 


Stars + WDs/BHs 


Total (Past SFR) 


Salpeter (1955) 


1.15 


1.30 


1.86 


Miller Scalo (1979) 


0.46 


0.60 


0.99 


Kennicutt (1983) 


0.46 


0.60 


1.06 


Scalo (1986) 


0.52 


0.61 


0.84 


Kroupa et al (1993) 


0.65 


0.76 


1.09 


Kroupa (2001) 


0.67 


0.83 


1.48 


Baldry & Glazebrook (2003) 


0.67 


0.86 


1.76 


Chabrier (2003) 


0.59 


0.75 


1.42 



masses of 0.1 and 120 Mq (Table 2). Stellar masses have been defined in 
various ways as represented by the 3 columns in Table 2. Typically we are 
interested in the observable stellar mass at a given time (i.e. main sequence 
and giant branch stars), but it is interesting to also compute tlic; total mass 
which is not in the interstellar medium, which includes that locked in evolved 
degenerate objects (white dwarfs and black holes). The most inclusive defini- 
tion of stellar mass (total) is the integral of the past star formation history. 
Depending on the definition, and chosen IMF, the uncertainties range almost 
over a factor of 4 for the most popularly- used functions, quite apart from the 
unsettling question of whether the form of the IMF might vary with epoch or 
type of object. 

Although the stellar mass function for a galaxy survey can be derived 
assuming a fixed mass/light ratio, the useful stellar density is that corrected 
for the fractional loss, R, of stellar material due to winds and supernovae. 
Only with this correction (if!=0.28 for a Salpeter IMF), does the present-day 
value represent the integral of the past star formation. 

Figure 18 shows the -fC-band luminosity and derived stellar mass function 
for galaxies in the 2dF redshift survey from the analysis of Cole et al (2001). 
if-band measures were obtained by correlation with the K <13.0 catalog 
obtained by the 2MASS survey. 

The integrated stellar density, corrected for stellar mass loss, is (Cole et 
al 2001): 

Qstars h = 0.0027 ± 0.00027 

for a Salpeter IMF, a value very similar to that derived independently by 
Fukugita & Peebles (2004). By comparison the local mass fraction in neutral 
HI -|- He I gas is: 

Qg^^ h = 0.00078 ± 0.00016 
Thus only 5% of all baryons are in stars with the bulk in ionized gas. 
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Fig. 18. (Left) Rest-frame iC-band luminosity function derived from the combi- 
nation of redshifts from the 2dF survey with photometry from 2MASS (Cole et al 
2001); Schechter parameter fits are shown. (Right) Derived stellar mass function 
assuming a Salpeter IMF corrected for lost material assuming _R=0.28 (see text). 

3.3 Diagnostics of Star Formation in Galaxies 

When significant redshift surveys became possible at intermediate and high 
redshift through the advent of multi-object spectrographs, so it became pos- 
sible to consider various probes of the star formation rate (SFR) at different 
epochs. As in the formalism for calculating the integrated luminosity density, 
PL, per comoving Mpc'^, so for a given population various diagnostics of on- 
going star formation can yield an equivalent global star formation rate psFR 
in units of M© yr"-'^ Mpc~^. 

Such integrated measures average over a whole host of important de- 
tails, such as differences in evolutionary behavior between luminous and sub- 
luminous galaxies and, of course, morphology. Moreover, in any survey at 
high redshift, only a portion of the population is rendered visible so uncertain 
corrections must be made to compare results at different epochs. The impor- 
tance of the cosmic star formation history, i.e. psfr{z), is it displays, in a 
simple manner, the epoch and duration of galaxy growth. By integrating the 
function, one should recover the present stellar density (§3.5). 

There are various probes of star formation in galaxies, each with its ad- 
vantages and drawbacks. Not only is there no single 'best' method to gauge 
the current star formation rate of a chosen galaxy, but as each probe samples 
the effect of young stars in different initial mass ranges, so each averages the 
star formation rate over a different time interval. If, as is often the case in the 
most energetic sources, the star formation is erratic or burst-like, one would 
not expect different diagnostics to give the same measure of the instantaneous 
SFR even for the same galaxies. 



Observations of the High Redshift Universe 



33 



Four diagnostics are in common use (see review by Kennicutt 1998). 

The rest-frame ultraviolet continuum (AA ~1250-1500A ) has the advan- 
tage of being directly connected to well- understood high mass (> SM©) 
main sequence stars. Large datasets are available for high redshift star- 
forming galaxies, including some to z ~6. Via the GALEX satellite and 
earlier balloon-borne experiments, local data is also available. The disad- 
vantage of this diagnostic lies in the uncertain (and significant) corrections 
necessary for dust extinction and a modest sensitivity to the assumed ini- 
tial mass function. Obscured populations are completely missed in UV 
samples. Kennicutt suggests the following calibration for the UV luminos- 
ity: 

SFR{MQyr-^) = lA10-'^^L^{ergs s'^ Hz-^) 

Nebular emission lines such as Ha and [O II] are also available for a range 
of redshifts {z < 2.5), for example as a natural by-product of faint redshift 
surveys. Gas clouds are photo-ionized by very massive (> IOMq) stars. 
Dust extinction can often be evaluated from higher order Balmer lines 
under various radiative assumptions depending on the escape fraction of 
ionizing photons. The sensitivity to the initial mass function is strong. 

SFR{MQyr-^) = 7.9 10-^2^(iJa) (ergs s"^) 

SFR{MQyr-^) = lA±0A10-'^^L{OII){ergs s''^) 

Far infrared emission (10-300 /J.m) arises from dust heated by young stars. 
It is clearly only a tracer in the most dusty systems and thus acts as a 
valuable complementary probe to the UV continuum. As we have seen in 
§2, luminous far infrared galaxies are also seen to high redshift. However, 
not all dust heating is due to young stars and the bolometric far infrared 
flux, Lpin, is needed for an accurate measurement. 

Radio emission, e.g. at 1.4GHz, is thought to arise from synchrotron emis- 
sion generated by relativistic electrons accelerated by supernova remnants 
following the rapid evolution of the most massive stars. Its great advan- 
tage is that it offers a dust-free measure of the recent SFR. Current radio 
surveys do not have the sensitivity to see emission beyond z ~1, so its 
promise has yet to b(^ fully explored. This process is also the least well- 
understood and calibrated. Sullivan et al (2001) discuss this point in some 
detail and conclude: 

SFR{MQyr-^) = l.llQ-'^^LiA{ergs s'^Hz-^) 
for bursts of duration >100 Myr. 
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The question of the time-dependent nature of the SFR is an important 
point (Sullivan et al 2000, 2001). For an instantaneous burst of star formation, 
Figure 19a shows the 'response' of the various diagnostics. Clearly if the SF 
is erratic on 0.01-0.1 Gyr timescales, each will provide a different sensitivity. 
Sullivan et al (2000) compared UV and Ha diagnostics for a large sample of 
nearby galaxies and found a scatter beyond that expected from the effects 
of dust extinction or observational error, presumably from this effect (Figure 
19b). 
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Fig. 19. Time dependence of various diagnostics of star formation in galaxies. (Left) 
sensitivities for a single burst of star formation. (Right) scatter in the UV and Ha 
galaxy luminosities in the local survey of Sullivan et al (2000); lines represent model 
predictions for various (constant) star formation rates and metallicities. It is claimed 
that some fraction of local galaxies must undergo erratic periods of star formation 
in order to account for the offset and scatter. 



In addition to the initial mass function (already discussed), a key uncer- 
tainty affecting the UV diagnostic is the selective dust extinction law. Over 
the wavelength range 0.3< A <1 /im, differences between laws deduced for 
the Milky Way, the Magellan clouds and local starburst galaxies (Calzetti 
et al 2000) are quite modest. Significant differences occur around the 2200 
A feature (dominant in the Milky Way but absence in Calzetti's formula) 
and shortward of 2000 A where the various formulae differ by ±2 mags in 
A{\)IE{B-V). 

3.4 Cosmic Star Formation - Observations 

Early compilations of the cosmic star formation history followed the field 
redshift surveys of Lilly et al (1996), Ellis et al (1996) and the abundance 
of U-band drop outs in the early deep HST data (Madau et al 1996). The 



Observations of the High Redshift Universe 



35 



pioneering papers in this regard include Lilly et al (1996), Fall et al (1996) 
and Madau et al (1996, 1998). 

Hopkins (2004) and Hopkins & Beacom (2006) have undertaken a valu- 
able recent compilation, standardizing all measures to the same initial mass 
function, cosmology and extinction law. They have also integrated the various 
luminosity functions for each diagnostic in a self-consistent manner (except 
at very high redshift). Accordingly, their articles give us a valuable summary 
of the state of the art. 

Figure 20 summarizes their findings. Although at first sight somewhat 
confusing, some clear trends are evident including a systematic increase in 
star formation rate per unit volume out to z ~1 which is close to (Hopkins 
2004): 

PSFRiz) oc (l + z)3-i 

A more elaborate formulate is fitted in Hopkins & Beacom (2006). 

There is a broad peak somewhere in the region 2< z <4 where the UV 
data is consistently an underestimate and the growing samples of sub-mm 
galaxies are valuable. The dispersion here is only a factor of ±2 or so, which 
is a considerable improvement on earlier work. We will return to the question 
of a possible decline in the cosmic SFR beyond z ~3-4 in later sections. 
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Fig. 20. Recent compilations of the cosmic star formation history. Circles are data 
from Hopkins (2004) color-coded by method: blue: UV, green: [O II], red: Ha//3, 
magenta: non-optical including sub-mm and radio. New data from Hopkins & Bea- 
com (2006), represented by various triangles, stars and squares, include Spitzer FIR 
measures (magenta triangles). The solid lines represent a range of the best fitting 
parametric form for z <1. 
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In their recent update, Hopkins & Beacom (2006) also parametrically fit 
the resulting Psfr{z) in two further redshift sections, beyond z and they 
use this to predict the growth of the absolute stellar mass density, p*, via 
integration (Figure 21). Concentrating, for now, on the reproduction of the 
present day mass density (Cole ct al 2001), the agreement is remarkably good. 

Although in detail the result depends on an assumed initial mass function 
and the vexing question of whether extinction might be luminosity-dependent, 
this is an important result in two respects: firstly, as an absolute comparison 
it confirms that most of the star formation necessary to explain the presently- 
observed stellar mass has already been detected through various complemen- 
tary surveys. Secondly, the study allows us to predict fairly precisely the epoch 
by which time half the present stellar mass was in place; this is zi = 2.0 ±0.2. 
In §4 we will discuss this conclusion further attempting to verily it by mea- 
suring stellar masses of distant galaxies directly. 




Fig. 21. Growth of stellar mass density, p*, with redshift obtained by direct inte- 
gration of a parametric fit to the cosmic star formation history deduced by Hopkins 

& Beacom (2006, sec Figure 20). The integration accurately reproduces the local 
stellar mass density observed by Cole et al (2001) and suggests half the present 
density was in place at ^:=2.0 ± 0.2. 



Observations of the High Redshift Universe 37 

3.5 Cosmic Star Formation - Theory 

As we have discussed, semi-analytical models have had a hard time repro- 
ducing and predicting the cosmic star formation history. Amusingly, as the 
data has improved, the models have largely done a 'catch-up' job (Baugh 
et al 1998, 2005a). To their credit, while many observers were still convinced 
galaxies formed the bulk of their stars in a narrow time interval (the 'primeval 
galaxy' hypothesis) , CDM theorists were the first to suggest the extended star 
formation histories now seen in Figure 20. 

A particular challenge seems to be that of reproducing the abundance 
of energetic sub-mm sources whose star formation rates exceed 100-200 Mq 
yr~^. Baugh et al (2005b) have suggested it may require a combination of qui- 
escent and burst modes of star formation, the former involving an initial mass 
function steepened towards high mass stars. Although there is much freedom 
in the semi-analytical models, recent models suggest zi ~ 1.3. By contrast, 
for the same cosmological models, hydrodynamical simulations (Nagamine et 
al 2004) predict much earlier star formation, consistent with zi ~ 2.0 — 2.5. 

The flexibility of these models is considerable so my personal view is that 
not much can be learned from these comparisons either way. It is more instruc- 
tive to compare galaxy masses at various epochs with theoretical predictions. 
Although we are still some ways from doing this in a manner that includes 
both baryonic and dark components, progress is already promising and will 
be reviewed in §4. 

3.6 Unifying the Various High Redshift Populations 

Integrating the various star-forming populations at high redshift to produce 
Figure 20 avoids the important question of the physical relevance and roles 
of the seemingly-diverse categories of high redshift galaxies. In the previous 
lecture (§2), I introduced three broad categories: the Lyman break (LBG), 
sub-mm and passively-evolving sources (DRGs) which co-exist over 1< z <3. 
What is the relationship between these objects? 

As the datasets on each has improved, we have secured important physical 
variables including masses, star formation rates and ages. We can thus begin 
to understand not only their relative contributions to the SFR at a given 
epoch, but the degree of overlap among the various populations. Several re- 
cent articles have begun to evaluate the connection between these various 
categories (Papovich et al 2006, Reddy et al 2005). 

A particularly valuable measure is the clustering scale, ro, for each popu- 
lation, as defined in §1.3. This is closely linked to the halo mass according to 
CDM and thus sets a marker for connecting populations observed at different 
epochs. Adelberger et al (1998) demonstrated the strong clustering, ro :ii3.8 
Mpc, of luminous LBGs at z ~3. Baugh et al (1998) claimed this was con- 
sistent with the progenitor halos of present-day massive ellipticals. The key 
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to the physical nature of LBGs depends the origin of their intense star for- 
mation. At z ~3, the bright end of the UV luminosity function is ~1.5 mags 
brighter than its local equivalent; the mean SFR is 45 Mq yr~^. Is this due 
to prolonged activity, consistent with the build up of the bulk of stars which 
reside in present-day massive ellipticals, or is it a temporary phase due to 
merger- induced starbursts (Somerville et al 2001). 

Shapley et al (2001, 2003) investigated the stellar population and stacked 
spectra of a large sample of z ~3 LBGs and find younger systems with intense 
SFRs are dustier with weaker Lya emission while outflows (or 'superwinds') 
are present in virtually all (Figure 22a). For the young LBGs, a brief period of 
elevated star formation seems to coincide with a large dust opacity hinting at 
a possible overlap with the sub-mm sources. During this rapid phase, gas and 
dust is depleted by outflows leading to eventually to a longer, more quiescent 
phase during which time the bulk of the stellar mass is assembled. 

If young dusty LBGs with SFRs ~ 3OOM0 yr~^ represent a transient 
phase, we might expect sub-mm sources to simply be a yet rarer, more extreme 
version of the same phenomenon. The key to testing this connection lies in 
the relative clustering scales of the two populations (Figure 22b). Blain et 
al (2004) find sub-mm galaxies are indeed more strongly clustered than the 
average LBGs, albeit with some uncertainty given the much smaller sample 
size. 
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Fig. 22. Connecting Lyman break and sub-mm sources. (Left) Correlation between 
the mean age (for a constant SFR) and reddening for a sample of z ~3 LBGs from 
the analysis of Shapley ot al (2001). The youngest LGBs seem to occupy a brief dusty 
phase limited eventually by the effect of powerful galaxy-scale outflows. (Right) The 
LBG-submm connection can be tested through reliable measures of their relative 
clustering scales, ro (Blain et al 2004). 
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Turning to the passively-evolving sources, although McCarthy (2004) pro- 
vides a valuable review of the territory, the observational situation is rapidly 
changing. For many years, CDM theorists predicted a fast decline with red- 
shift in the abundance of red, quiescent sources. Using a large sample of 
photometrically-selected sources in the COMBO-17 survey. Bell ct al (2004) 
claimed to see this decline in abundance by witnessing a near-constant lumi- 
nosity density in red sources to ^; ~1 (Figure 23a). The key point to under- 
stand here is that a passively-evolving galaxy fades in luminosity so that the 
red luminosity density should increase with redshift unless the population is 
growing. Bell et al surmises the abundance of red galaxies was 3 times less at 
z cil as predicted in early semi-analytical models (Kauffmann et al 1996). 

By contrast, the Gemini Deep Deep Survey (Glazebrook et al 2004) finds 
numerous examples of massive red galaxies with ^; > 1 in seeming contradic- 
tion with the decline predicted by CDM supported by Bell et al (2004). Of 
particular significance is the detailed spectroscopic analysis of 20 red galaxies 
with z ~1.5 (McCarthy et al 2004) whose inferred ages are 1.2-2.3 Gyr imply- 
ing most massive red galaxies formed at least as early as z ~2.5-3 with SFRs 
of order 300-500 Mq yr~^. Could the most massive red galaxies at z ~1.5 then 
be the descendants of the sub-mm population? One caveat is that not all the 
stars whose ages have been determined by McCarthy et al need necessarily 
have resided in single galaxies at earlier times. The key question relates to the 
reliability of the abundance of early massive red systems. Using a new color- 
selection technique, Kong et al (2006) suggest the space density of quiescent 
systems with stellar mass > IO^^Mq at z ~1.5-2 is only 20% of its present 
value. 

As we will see in §4, the key to resolving the apparent discrepancy between 
the declining red luminosity density of Bell et al and the presence of massive 
red galaxies at z ~1-1.5, lies in the mass-dependence of stellar assembly (Treu 
et al 2005). 

Finally, a new color-selection has been proposed to uniformly select all 
galaxies lying in the strategically-interesting redshift range 1.4< z <2.5. 
Daddi et al (2004) have proposed the 'BzK' technique, combining (z — K) 
and [B — z) to locate both star-forming and passive galaxies with z >1.4; 
such systems are termed 'sBzK' and 'pBzK' galaxies respectively. Reddy et 
al (2005) claim there is little distinction between the star-forming sBzK and 
Lyman break galaxies - both contribute similarly to the star formation density 
over 1.4< z <2.6 and the overlap fractions are at least 60-80%. 

More interestingly, both Reddy et al (2005) and Kong et al (2006) suggest 
significant overlap between the passive and actively star-forming populations. 
Kong et al find the angular clustering is similar and Reddy et al find the 
stellar mass distributions overlap. 
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Fig. 23. Left: The rest-frame blue luminosity density of 'red sequence' galaxies 
as a function of redshift from the COMBO-17 analysis of Bell et al (2004). Since 
such systems should brighten in the past the near-constancy of this density implies 
3 times fewer red systems exist at z ~1 than locally - as expected by standard 
CDM models. Right: Composite spectra of red galaxies with 1.3< z <1.4 (blue), 
1.6< z <1.9 (black) from the Gemini Deep Deep Survey (McCarthy et al 2004). 
A stellar synthesis model with age 2 Gyr is overlaid in red. The analysis suggests 
the most massive red galaxies with z ~5 formed with spectacular SFRs at redshifts 
z ~2.4-3.3. 

3.7 Lecture Summary 

Clearly multi- wavelength data is leading to a revolution in tracking the history 
of star formation in the Universe. Because of the vagaries of the stellar initial 
mass function, dust extinction and selection biases, we need multiple probes 
of star formation in galaxies. 

The result of the labors of many groups is a good understanding of the 
comoving density of star formation since a redshift z ~3. Surprisingly, the 
trends observed can account with reasonable precision for the stellar mass 
density observed today. The implication of this result is that half the stars we 
see today were in place by a redshift z ~2. 

What, then, are we to make of the diversity of galaxies we observe during 
the redshift range {z ~2-3) of maximum growth? Through detailed studies 
some connections are now being made between both UV-emitting Lyman 
break galaxies and dust-ridden sub- mm sources. 

More confusion reigns in understanding the role and decline with redshift 
in the contribution of passively-evolving red galaxies. Some observers claim 
a dramatic decline in their abundance whereas others demonstrate clear evi- 
dence for the presence of a significant population of old, massive galaxies at 
z ~1.5. We will return to this enigma in §4. 
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4 Stellar Mass Assembly 

4.1 Motivation 

Although wc may be able to account for the present stellar mass density by 
integrating the comoving star formation history (§3), this represents only a 
small step towards understanding the history of galaxy assembly. A major 
limitation is the fact that the star formation density averages over a range of 
physical situations and luminosities; we are missing a whole load of important 
physics. As we have seen, a single value for the stellar mass density, e.g. 
p^{z = 2), is useful when considered as a global quantity (e.g. compared to an 
equivalent estimate at z=0), but it does not describe whether the observed 
star formation is steady or burst-like in nature, or even whether the bulk of 
activity within a given volume arises from a large number of feeble sources 
or a small number of intense objects. Such details matter if we are trying to 
construct a clear picture of how galaxies assemble. 

Of course we could extend our study of time-dependent star formation 
to determine the distribution functions of star formation at various epochs 
(e.g. the UV continuum. Ha or sub- mm luminosity functions), but making 
the integration check only at z=0 is second-best to measuring the assembled 
mass and its distribution function at various redshifts. This would allow us to 
directly witness the growth rate of galaxies of various masses at various times 
and, in some sense, is a more profound measurement, closer to theoretical 
predictions. 

Ideally we would like to measure both baryonic and non-baryonic masses 

for large numbers of galaxies but, at present at least, we can only make dynam- 
ical or lensing-based total mass estimates for specific types of distant galaxy 
and crude estimates for the gaseous component. The bulk of the progress made 
in the last few years has followed attempts to measure stellar masses for large 
populations of galaxies. We will review their achievements in this lecture. 

4.2 Methods for Estimating Galaxy Masses 

What are the options available for estimating galaxy masses of any kind at 
intermediate to high redshift? Basically, we can think of three useful methods. 

• Dynamical methods based on resolved rotation curves for recognisable disk 
systems (Vogt et al 1996,1997) or stellar velocity dispersions for pressure- 
supported spheroidals (van Dokkum & Ellis 2003, Treu et al 2005). These 
methods only apply for systems known or assumed a priori to have a 
particular form of velocity field. Interesting constraints are now available 
for a few hundred galaxies in the field and in clusters out to redshifts of 
z ~1.3. Key issues relate to biases associated with preferential selection of 
systems with 'regular' appearance and how to interpret mass dynamically- 
derived over a limited physical scale (c.f. Conselice et al 2005). In the 
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absence of resolved rotation curves, sometimes emission linewidths are 
considered a satisfactory proxy (Newman & Davis 2000). 

• Gravitational lensing offers the cleanest probe of the total mass distribu- 
tion but, as a geometric method, is restricted in its application to compact, 
dense lenses (basically, sphcroidals) occupying cosmic volumes typically 
half way to those probed by faint star-forming background field galaxies. 
In practice this means ziens <1- Even so, by sifting through spectra of 
luminous red galaxies in the SDSS survey and locating cases where an 
emission line from a background lensed galaxy enters the spectroscopic 
fiber, Bolton et al (2006) have identified a new and large sample of Ein- 
stein rings enabling us to gain valuable insight into the relative distribution 
of dark and visible mass over 0< z <1. 

• Stellar masses derived from near-infrared photometry represents the most 
popular technique in use at the current time. The idea has its origins in the 
recognition (Broadhurst, Ellis & Glazebrook 1992, Kauffmann & Chariot 
1998) that the rest-frame if-band luminosity of a galaxy is less affected by 
recent star formation than its optical equivalent (Figure 16), and thus can 
act as a closer proxy to the well-established stellar population. A procedure 
for fitting the rest-frame optical-infrared spectral energy distribution of a 
distant galaxy, deriving a stellar mass/light ratio (M/Lk) and hence the 
stellar mass if the redshift is known, was introduced by Brinchmann & 
Ellis (2000). The popularity of the technique follows from the fact it can 
easily be applied to large catalogs of galaxies in panoramic imaging surveys 
and extended to very high redshift if IRAC photometry is available. The 
main difficulties relate to the poor precision of the method, particularly if 
the same photometric data is being used to estimate the redshift (Bundy 
et al 2006), plus degeneracies arising from poor knowledge of the past star 
formation history (Shapley et al 2005, van der Wei et al 2006). 

For many galaxies, an important and usually neglected component is the 
mass locked up in both ionized and cool gas. In nearby systems amenable to 
study of hot ionized gas (from nebular emission lines) and its usually dom- 
inant cooler neutral component (probed by 21cm studies), as much as 20% 
of the baryonic mass of a luminous star-forming galaxy can be found in this 
form (Zwaan et al 2003). At present, it is not possible to routinely use radio 
techniques to reliably estimate gaseous masses of distant galaxies although 
approximate gas masses have been derived assuming the projected surface; 
density of nebular emission correlates with the gas mass within some mea- 
sured physical scale (Erb et al 2006). 

4.3 Results: Regular Galaxies 0< z <1.5 

Because of the simplicity of their stellar populations, velocity fields and the 
lack of confusing gaseous components, rather more is known about the mass 
assembly history of ellipticals than for spirals. Concerning ellipticals, one of 
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the key challenges is separating the age of the stars from the age of the as- 
sembled mass. 

The Fundamental Plane (FP, Dressier et al 1987, Djorgovski & Davis 1987, 
Bender, Burstein & Faber 1992, Jorgensen et al 1996) represents an empirical 
correlation between the dynamical mass (via the central stellar velocity dis- 
persion (To), the effective radius (Re) and light distribution (via the enclosed 
surface brightness he) for ellipticals, viz: 

logRE — alogaQ + b ^e + c 

For example, with a in km sec~^. Re in kpc and /^b in mags arcsec"^ in 
the B band, a=1.25, 6=0.32, c=-8.970 for h=0.65. 

These observables define an effective dynamical mass Me oc a^RE/G 
which correlates closely with that deduced from lensing (Treu et al 2006). 
Deviations from the local FP at a given redshift z can be used to deduce the 
change in mass/light ratio Alog{M / L). 

The most comprehensive studies of the evolving FP come from two inde- 
pendent and consistent studies of field spheroidals to 2 ~1 (Treu et al 2005, 
van der Wei et al 2005). The evolution in mass/light ratio Alog{M/L) de- 
duced from the GOODS-N survey of Treu et al (2005) is shown in Figure 
24. These authors find as little as 1-3% by stellar mass of the present-day 
population in massive (> lO^^'^M©) galaxies formed since z=1.2, whereas for 
low mass systems (< IQ^^Mq) the growth fraction is 20-40%. This result, 
confirmed independently by van der Wei et al (2005), is an important illustra- 
tion of the mass-dependent growth in galaxies with the most massive systems 
shutting off earliest. 

Of course, one should not confuse the age of stars, as probed by the FP, 
with the age of the assembled mass, van Dokkum (2005) has argued that 
if spheroidals preferentially merge with similar gas-poor systems (a process 
called 'dry mergers') the FP analyses could well indicate early eras of major 
star formation even though the bulk of the assembled mass in individual sys- 
tems occurred at z <l|3. Bell et al (2006) and Tran et al (2005) have cataloged 
individual cases of dry mergers in both field and cluster samples, respectively. 
Their occurrence is not in dispute; however, only via morphological or other 
measures of the global mass assembly can the role of dry mergers as a major 
feature of galaxy assembly be addressed. 

Related insight into this problem arises from the relative distributions of 
baryonic and dark matter deduced from the combination of lensing and stel- 
lar dynamics for the recently-located SDSS-selected Einstein rings (Treu et 
al 2006, Figure 25). Although it might be thought that lensing preferentially 
selects the most massive and compact sources, Treu et al compare the FP of 

* The reason observers have gone somewhat out of their way to consider such 
complicated scenarios is because late assembly of massive spheroidals was, until 
recently, a fundamental tenet of the CDM hierarchy to be salvaged at all costs 
(KaufFmann et al 1996). 
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Fig. 24. Change in B-band mass/light ratio with redshift deduced from the Keck 
dynamical survey of over 150 field spheroidals in GOODS-N (Treu et al 2005). The 
plot shows the change in mass/light ratio for galaxies of different effective dynamical 
masses (red: massive, cyan: intermediate, blue: low mass). Curves illustrate the 
change in mass/light ratio expected as a function of redshift for galaxies assembled 
monolithically with simply evolving stellar populations since a redshift of formation 
Zf. Clearly the stars in the most massive galaxies largely formed at high redshift 
whereas assembly continued apace in lower mass systems. 

such lenses with those in the larger field sample (Fig. 25) and deduce other- 
wise. An important result from the study of the first set of such remarkable 
lenses is how well the total mass profile can be represented by an isothermal 
form with mass tracing light, viz: 

r _ 
Ptot = P( — ) 
ro 

Even over 0< z <1, the mass slope 7 is constant at 2.0 to within 2% 
precision indicating rather precise coUisional coupling of dark matter and gas. 

Sadly, far less is known about the mass assembly history of regular spirals. 
The disk scaling law equivalent to the FP for ellipticals, the TuUy-Fisher 
relation (TuUy & Fisher 1977) which links rotational velocity to luminosity 
gives ambiguous information without additional input. Modest evolution in 
the TF relationship was deduced from the pioneering Keck study of z ~1 
spirals by Vogt et al (1997) but this could amount to ~0.6 mag of B-band 
luminosity brightening in sources of a fixed rotational velocity to z ~1, or 
more enhancement if masses were reduced. 

Additional variables capable of breaking the degeneracy between dynam- 
ical mass and luminosity include physical size and stellar mass. Lilly et al 
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Fig. 25. Distribution of the gradient gamma in the total mass profile with redshift 
as deduced from the combination of stellar dynamics and lensing for ellipticals in the 
SLAG survey (Treu et al, in prepn). The remarkable constancy of the profile slope 
indicates massive relaxed systems were in place at z ~1 and that 'dry merging' 
cannot be a prominent feature of the assembly history of large galaxies. 



(1998) examined the size-luminosity Telaiion for several hundred disks to z ~1 
in a HST-classified redshift survey sample (see Sargent et al 2006 for an up- 
date) and found no significant growth for the largest systems. Conselice et al 
(2005) correlated stellar and dynamical masses for ~100 spirals with resolved 
dynamics in the context of a simple halo formulation. Although their deduced 
halo masses must be highly uncertain, they likewise deduced that growth must 
be modest since z ~1, occurring in a self-similar fashion for the baryonic and 
dark components. 

4.4 Stellar Masses from Multi-Color Photometry 

Bundy et al (2005, 2006) give a good summary and critical analysis of the now 
well-established practice of estimating stellar masses from multi-color optical- 
infrared photometry. Figure 26 gives a practical illustration of the technique 
where it can be seen that even for low z galaxies with good photometry, 
the precision in mass is only ~ ±0.1-0.2 dex. In most cases, even random 
uncertainties are at the ±0.2-0.3 dex level and systematic errors are likely to 
be much higher. 

Since analysing stellar mass functions is now a major industry in the com- 
munity, it is worth spending some time considering the possible pitfalls. A 
significant fraction of the large datasets being used are purely photometric, 
with both redshifts and stellar masses being simultaneously deduced from 
multi-color photometry (Drory et al 2002, Bell et al 2004, Fontana et al 2004). 
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Fig. 26. Deriving stellar masses for galaxies of known redshift via multi-color optical 
photometry (after Bundy et al 2006). (Left) Rest-frame spectral energy distribution 

of two DEEP2 galaxies of known redshift from broad-band BRIK photometry with 
a fit deduced by fitting from a stellar synthesis library. The fit yields the most likely 
mass/light ratio of the observed population. (Right) Likelihood distribution of the 
stellar mass for the same two galaxies derived by Bayesian analysis. Dashed lines 
indicate the eventual range of permitted solutions. 

Few large surveys have extensive spectroscopy from which to check that this 
procedure works. 

Bundy et al (2006) address this important question in the context of 
their extensive DEEP2 spectroscopic sample by artificially 'switching off' the 
knowledge of the spectroscopic z: how does the stellar mass deduced when si- 
multaneously deriving the photometric redshift from the same data compare 
with that derived when the spectroscopic value is externally input? Figure 27 
illustrates that a significant component of error is one beyond that expected 
solely from the error in the derived redshift (as measured from the Zgpec vs. 
Zphoto comparison). 

A second restriction in many stellar mass determinations is the absence 
of any near-infrared photometry. This may seem surprising given the classic 
papers (Kauffmann & Chariot 1998, Brinchmann & Ellis 2000) stressed its 
key role. However, panoramic near-infrared imaging is much more expensive 
in telescope time as few observatories, until recently, had large format infrared 
cameras. Inevitably, some groups have attempted to get by without the near- 
infrared data. 
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Fig. 27. Error in log steUar mass when the known spectroscopic redshift is ignored 
(Bundy et al 2005). The shaded region defines the expected scatter in mass arising 
solely from that in the photometric-spectrosopic redshift comparison. Clearly in 
addition to this component, some catastrophic failures are evident. 



Figure 28 illustrates how the precision degrades when the K band data is 
dropped in the stellar mass fit (Bundy et al 2006 jl. Although the systematic 
error is contained, the noise increases significantly for redshifts z >0.7 where 
the optical photometry fails to adequately sample the older, lower mass stars. 

Shapley et al (2005) explore a further uncertainty, which is particularly 
germane to the analysis of their z ~2 sources. Using IRAC to represent the 
rest- frame infrared at these redshifts, they consider the role that recent bursts 
of star formation might have on the deduced stellar masses. Although bursts 
predominantly affect short-wavelength luminosities, one might imagine little 
effect at the longer wavelengths (c.f. Kauffmann & Chariot 1998). Shapley et 
al consider a wider range of star formation histories and show that the derived 
stellar mass depends less strongly on the long wavelength luminosity L (4.5^m) 
than expected. The scatter is consistent with variations in mass/light ratios 
of xl5 (Figure 29). The strong correlation between mass and optical- infrared 
color emphasizes the importance of secondary activity. 



* A similar analysis was conducted by Shapley et al (2005) at z ~2, excluding the 
relevant IRAC data. 
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Fig. 28. Effect of deriving stellar masses from optical photometry alone (Bundy 
et al 2006). (Left) Difference in masses derived using BRIK and BRI photometry. 
(Center) Distribution of uncertainties deduced from Bayesian fitting for the BRI 
fits. (Right) as center for the BRIK fitting. Although the systematic offset is small, 
the uncertainty in deduced stellar mass increases dramatically for galaxies with 
z >0.7. 
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Fig. 29. Left: Stellar masses versus observed IRAC 4.5pim luminosities for the 
sample of z ~2 galaxies studied by Shapley et al (2005). Allowing a wide range of 
star formation histories illustrates a weaker dependence of mass on long wavelength 
luminosity than seen in earlier tests at lower redshift based on simpler star formation 
histories. The range in rest-frame 1.6/im mass/light ratio is ~ xl5. (Right) The 
existence of a fairly tight correlation between mass and R — K color demonstrates 
the importance of considering secondary bursts of star formation in deriving stellar 
masses. 
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4.5 Results: Stellar Mass Functions 0< z <1.5 

Noting the above uncertainties, we now turn to results compiled from the 
stellar mass distributions of galaxies undertaken in both spectroscopic and 
photometric surveys. A recent discussion of the various results can be found 
in Bundy et al (2006). 

Brinchmann & Ellis (2000) were the first to address the global evolution of 
stellar mass in this redshift range using a morphological sample of 350 galaxies 
with spectroscopic redshifts and infrared photometry. The small sample did 
not permit consideration of detailed stellar mass functions, but the integrated 
mass density was partitioned in 3 redshift bins for spheroidals, spirals and 
irregulars. Surprisingly little growth was seen in the overall mass density from 
^; ~1 to today; the strongest evolutionary signal seen is a redistribution of 
mass amongst the morphologies dominated by a declining mass density in 
irregulars with time. 

The ratio of the galactic stellar mass Mgtar to the current star formation 
rate SFR, e.g. as deduced spectroscopically or from the rest-frame optical 
colors, is sometimes termed the specific star formation rate, R. This quantity 
allows us to address the question of whether galaxies have been forming stars 
for a significant fraction of the Hubble time, at a rate commensurate with 
explaining their assembled mass. A low value for R implies a quiescent object 
whose growth has largely ended; the mean stellar age is quite largo. A high 
value of i? implies an active object which has assembled recently. A frequently- 
used alternative is the 'doubling time' - that period over which, at the current 
SFR, the observed stellar mass would double. This time would be quite short 
for active objects. 

Cowie et al (1996), Brinchmann & Ellis (2000), and most recently Juneau 
et al (2005), found a surprising trend whereby most massive galaxies over 
z ~0.5-1.5 are quiescent, having presumably formed their stars well before 
z ~2, whereas low mass galaxies remain surprisingly active. The term 'down- 
sizing' - a signature of continued growth in lower mass galaxies after that in 
the high mass galaxies has been completed - was first coined by Cowie et 
al (1996) and has been used rather loosely in the recent literature to imply 
any signature of anti- hierarchical activity. In particular, it is important to 
distinguish between downsizing in star formation activity, which prcsmnably 
represents some physical process that permits continued star formation in 
lower mass systems when that in massive galaxies has conchided, from down- 
sizing in mass assembly, a truly 'anti-hierarchical' process whereby new mass 
is being added to lower mass galaxies at later times (see discussion in Bundy 
et al 2006). 

Before trying to understanding in more detail what causes this mass- 
dependent star formation, it is worth returning to the issue of dry merg- 
ers raised in §4.3. A 'downsizing' signature was also seen in the growth of 
spheroidal galaxies as analysed by their location on the evolving Fundamen- 
tal Plane (Treu et al 2005, van der Wei et al 2005). However, van Dokkum 
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(2005) argued that while this may reflect older stars in the more massive 
galaxies, the assembled mass may still be young if merging preferentially oc- 
curs between quiescent objects. The only way to test this hypothesis is to 
directly measure the growth rate in spheroidal systems. 

Using the deeper and more extensive sample of galaxies available in the 
GOODS fields, Bundy et al (2005) produced type-dependent stellar mass 
functions (Figure 30). Here, for the first time, one can see the morphological 
evolution detected by Brinchmann & Ellis largely arises via the transforma- 
tion of intermediate mass (~5. 10^" - 2.10^^M©) irregulars and spirals into 
spheroidals. If merging is responsible for this transformation, it is predomi- 
nantly occurring between gas-rich and active systems even at the very highest 
masses. The bulk of the evolution below z ~1 is simply a redistribution of 
star formation activity, perhaps as a result of mergers or feedback processes. 

Interestingly, in Figure 30 there is almost no change in the total mass 
function with time above 5. IO^^Mq, suggesting little growth in the mass 
spheroidals that dominate the high mass end. However, it worth remembering 
that the two GOODS fields are limited in size (0.1 deg^) and suffer from 
cosmic variance at the 20% level at high redshift increasing to 60% at low 
redshift (Bundy et al 2005). 

Stellar mass functions for a much larger sample of field galaxies of known 
redshift have been analyzed by Bundy et al (2006) utilizing the combination 
of extensive spectroscopy and Palomar i^-band imaging in four DEEP2 fields 
(totalling 1.5 deg^). This sample has the benefit of being much less affected 
by cosmic variance although, as there is no complete coverage with HST, mor- 
phological classifications are not possible. Color bimodality has been analyzed 
in the DEEP2 sample (Weiner et al 2005) and Bimdy et al use the rest-frame 
U — B color and spectroscopic [O II] equivalent width to separate quiescent 
and active galaxies. 

Figure 31 shows a direct comparison of the integrated stellar mass func- 
tions from this large survey alongside other, less extensive surveys, most of 
which are based only on photometric redshifts. Although each is variance lim- 
ited in different ways, one is struck again by the quite modest changes in the 
abundance of massive galaxies since z ~1. 

Next we consider (Figure 32) the stellar mass functions for the quiescent 
and active star-forming galaxies independently, partitioned according to the 
rest-frame U—B color. The surprising result here is the existence of a threshold 
or quenching mass, Mq above which there are no active systems. This is 
implied independently in Figure 30 where, within the redshift bin 0.55 < z < 
0.8, there are no spirals or irregulars with mass > 10^^ Mq. Interestingly, there 
is also a modest downward transition in Mq with time. 

A remarkably consistent picture emerges from these studies. Over the red- 
shift range 0< z <1.5, the stellar mass function has changed very little at 
the high mass end. Substantial growth, in terms of a mass doubling on ~5-8 
Gyr timescalc, is only possible for galaxies with stellar masses below lO"'^" Mq. 
Above this mass, the basic evolutionary signal is a quenching of star forma- 




Fig. 30. Stellar mass function arranged by morphology in two redshift bins from the 
analysis of Bundy et al (2005) for h—1. The solid curve in both plots represents the 
present day mass function from the 2dF survey (Cole et al 2003). Type-dependent 
mass functions are color-coded with black representing the total. The solid line 
connecting the filled black circles represents the sample with spectroscopic redshifts, 
the dotted line connecting the open black diamonds includes masses derived from 
sources with photometric redshifts. The vertical dashed lines represent completeness 
limits for all types. 



tion in well-established systems. This quenching progresses toward lower mass 
systems at later times, consistent with the mass-dependent trends seen in the 
ages of stellar populations in the quiescent spheroidals. 

The physical origin of the quenching of star formation, the fundamental 
origin of the various downsizing signals, is unclear. Although the merger- 
induced production of active galactic nuclei may lead to the temporary ex- 
pulsion and heating of the gaseous halos that surround galaxies (Springel et 
al 2005, Croton et al 2006), key tests of this hypothesis include sustaining the 
quenching, the weak environmental dependence of the observational trends 
and the surprisingly clear redshift dependence of the effect (Figure 32). 
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Fig. 31. The evolving stellar mass function from the comparison of Bundy et al 
(2006). 



4.6 Results: Stellar Mass Functions z>1.5 

With current facilities, the stellar mass data beyond z ~1.5 generally probes 
only the higher mass end of the distribution and relies on photometric rather 
than spectroscopic data. Nonetheless, the results emerging have received as 
much attention as those at lower redshift. Unlike the complexities of under- 
standing downsizing and the redistribution of mass and morphology in the 
z <1 data, the basic question at stake here is simply whether the abundance 
of massive z ^2 galaxies is larger than expected in the standard model. 

Testing the decline with redshift in the comoving abundance of, say, sys- 
tems with stellar mass greater than 10^^ Mq, expected in A CDM has been a 
frustrating story for the observers for two reasons. Firstly, the most massive 
systems are rare and clustered, and so determining reliable density estimates 
beyond ^; ~1.5 has required panoramic deep infrared data which has only re- 
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Fig. 32. The evolving stellar mass function split into quiescent (red) and active 
(blue) sources from the survey of Bundy et al (2006). Uncertainties arising from 
counting statistics and errors in the stellar masses are indicated by the associated 
shading. The thin solid line represents the local 2dF function (Cole et al 2003) and 
the dark and light grey regions represent total mass functions using only spectro- 
scopic and including photometric redshifts respectively. 

cently arrived. For a review of the early observational efforts see Benson, Ellis 
& Menanteau (2002) 

Secondly, there has been considerable confusion in the theoretical litera- 
ture on the expected rate of decline in > 10^"'^ Mq systems. Early predictions 
(Kauffmann, Chariot & White 1996) claimed a 3-fold decline in the comoving 
abundance to z ~1 in apparent agreement with the large photometric sample 
analyzed by Bell et al (2004). However, careful comparisons of independent 
semi-analytical predictions (Kauffmann et al 1999, Cole et al 2000) reveal 
substantial differences (by an order of magnitude for the same world model) 
in the rate of decline even to z ~1 (Benson, Ellis & Menanteau 2002). In re- 
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ality, the predictions depend on many parameters where differences can have 
a large effect for the region of the mass function where the slope dN / dM is 
steep. 

Dickinson et al (2003) undertook a pioneering study to evaluate the growth 
of stellar mass with photometric redshift using deep H <26.5 NICMOS data 
in the Hubble Deep Field North. Although a tiny field (5 arcmin^), the work 
was the first to demonstrate the existence of 10^^ Mq galaxies at z ~1.5-2 
as well as the challenges of estimating reliable abundances. A similar HDF-S 
analysis was undertaken by Rudnick et al (2003). 

A more representative area of 4x30 arcmin was probed spectroscopically 
in the Gemini Deep Deep Survey (Glazcbrook et al 2004), albeit to a much 
shallower depth. These authors claimed a 'surprising' abundance of 10^^ Mq 
systems to z ~2. Although an absolute comparison of the abundances with 
semi-analytical predictions is not likely to be illuminating for the reasons 
mentioned above, the redshift-dependent mass growth rate over 1< z <2, 
derived empirically from both the stellar masses and the integrated star for- 
mation rate, seems much slower than expected in the semi-analytical models 
(Granato et al 2000). Of particular interest is the fact that the mass growth 
rate seems not depend strongly on the mass range in question in apparent 
disagreement with the model predictions (Figure 33). 

Although many of the observers associated with the Gemini Deep Deep 
Survey have claimed the abundance and slow growth rate of massive galaxies 
over 1.5< z <2 poses a crisis for the standard model, alongside the puzzling 
'anti-hierarchical' behavior observed for z <1.5, the plain fact is that there is 
considerable uncertainty in the semi-analytical predictions. 

4.7 Quiescent Galaxies with 2< z <3 

Finally, it is illustrative to consider the dramatic effect that infrared data 
from panoramic ground-based cameras and the Spitzer Space Telescope is 
having, not only on our knowledge of the distribution of stellar masses at 
high redshift, but also how stellar mass is distributed among quiescent and 
star-forming populations. Until recently, there was widespread belief that the 
bulk of the star formation in this era, and probably a significant fraction of the 
stellar mass, lay in the Lyman break population. By contrast, van Dokkum 
et al (2006) have examined the rest-frame U — V colors of a sample of 300 
> 10^-^ Mq galaxies with 2< Zptoto < 3 (Figure 34) and claim almost 80% 
are quiescent. The definition of 'quiescent' here is somewhat important to 
get correct, given it now emerges that many of the originally-selected distant 
red galaxies with J — K >2.3 (§2.6) turn out to have quite respectable star 
formation rates. Infrared spectroscopy sensitive to Ha emission and the 4000A 
break is making great strides in clarifying this question (Kriek et al 2006). 




Fig. 33. The mass assembly history derived from stellar masses in the Gemini 
Deep Deep spectroscopic survey (Glazebrook et al 2004). Colored data points refer 
to observed cumulative densities p ( > M) for various mass ranges and the green 
shaded region to estimates derived from the integrated star formation history. The 
colored lines refer to predicted growth rates according to the GALFORM semi- 
analytic models. 
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Fig. 34. The distribution of rest-frame U-V color for a Spitzer-selected sample of 
2 < z <3 galaxies with stellar masses > 10^^ Mq from the analysis of van Dokkum 
et al (2006). Almost 80% of this sample in a quiescent state consistent with having 
completed the bulk of their assembly. 
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4.8 Lecture Summary 

In this lecture we have demonstrated important new techniques for estimating 
the stellar masses of all types of galaxies to impressive high redshift {z ~3). 
The reliability of these techniques is improved greatly by having spectro- 
scopic, rather than photometric, redshifts and, for z > 0.7, by the addition 
of infrared photometry. These techniques augment and extend more precise 
measures available for restricted classes of galaxies - such as the Fundamental 
Plane for pressure-supported spheroidals, and the TuUy-Fisher relationship 
for rotationally-supported disks. 

These various probes point to a self-consistent, but puzzling, description 
of mass assembly since z ~1.5. During this era, massive galaxies have hardly 
grown at all - most of the star formation and rapid growth is occurring in 
systems whose masses arc less than 10^" Mq. Remarkably, the star formation 
appears to be quenched above a certain threshold mass whose value, in turn, 
is declining with time. 

Energetic sources such as supernovac or active nuclei may be responsi- 
ble for this 'downsizing' signature but further work is needed to verify both 
the weak environmental trends seen in the observations and the redshift- 
dependent trends in the threshold mass. 

Beyond z ~1.5, the number and distribution of massive (> lO^"*^ Mq galax- 
ies has led to some surprises. Although the sheer abundance may not be a 
problem for contemporary models, the fact that so many have apparently com- 
pleted their star formation is more challenging and consistent with the slow 
growth rate observed at later times. The observational situation is rapidly 
developing but consistent with the presence of a surprisingly abundant and 
mature population of massive galaxies by z ~3. 
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5 Witnessing the End of Cosmic Reionization 
5.1 Introduction - Some Weighty Questions 

Wc now turn to the exciting and rapidly developing area of imderstanding 
cosmic reionization. This event, which marks the end of the so-called 'Dark 
Ages' when the intergalactic medium became transparent to ultraviolet pho- 
tons, was a landmark in cosmic history. In some ways the event might be 
considered as important as the epoch of recombination which isolates the for- 
mation of the hydrogen atom, or the asssociated surface of last scattering 
when photons and baryonic matter decouple. In the case of the era of cosmic 
reionization, although we cannot yet be sure, many believe we are isolating 
that period when the first sources had sufficient output to contribute to the 
energy balance of the intergalactic medium. Even though some early lumi- 
nous forerunners might be present, the epoch of reionization can be directly 
connected with cosmic dawn for starlight. 

It seems an impossible task to give an authoritative observational account 
of how to probe this era. So many issues are complete imponderables! When 
did reionization occur? Was it a gradual event made possible by a complex 
time sequence of sources, or was there a spectacular synchronized moment? 
Can we conceive of an initial event, followed by recombination and a second 
phase? 

What were the sources responsible? History has shown the naivety of as- 
tronomers in assuming a single population to be responsible for various phe- 
nomena - usually some complex combination is the answer. And, perhaps most 
ambitiously, what is the precise process by which photons escape the sources 
and create ionized regions? 

Four independent observational methods are helpful in constraining the 
redshift range where we might search for answers to the above questions. In 
this lecture, we will explore how these work and the current (and rapidly 
changing) constraints they offer. These are: 

• Evolution in the optical depth of Lyman series absorption lines in high 

redshift quasars (Fan ct al 2006a, b). Although QSOs are only found to 
redshifts of z ~6.5, high quality spectroscopic data gives us a glimpse of 
a potential change in the degree of neutrality of the IGM beyond z ~5.5 
such as might be expected if reionization was just ending at that epoch. 

• The ubiquity of metal absorption lines in the IGM as probed again by 
various spectra of the highest redshift QSOs. Carbon, in particular, is 
only produced in stellar nuclei and its presence in all sightlines to z ~5- 
6 (Songaila 2004) is highly indicative of a widespread process of earlier 
enrichment from the first generation of supernovae. 

• The large angular scale power in the temperature - polarization cross- 
correlation function seen in the microwave background (Kogut et al 2003, 
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Spergel et al 2006). This signal is produced by electron scattering fore- 
ground to the CMB. Depending on models of large scale structure, the op- 
tical depth of scattering gives some indication of the redshift range where 
the ionized particles lie. 
• The stellar mass density in assembled sources at rcdshifts z ~5-6, about 1 
billion years after the Big Bang, as probed by the remarkable combination 
of HST and Spitzer (Stark fc Ellis 2005, Stark et al 2006a, Yan et al 2006, 
Eyles et al 2006). A census of the mass in stars must be the integral of 
past activity which can be compared with that necessary to reionize the 
Universe. 

5.2 The Gunn-Peterson Test and SDSS QSOs 

In a remarkable paper, long before QSOs were located at redshifts beyond 
2.5, James Gunn and Bruce Peterson realized (Gunn & Peterson 1965) that 
the absence of broad troughs of hydrogen absorption in the spectra of QSOs 
must indicate intergalactic hydrogen is ionized. They postulated a future test 
whereby the spectra of QSOs of successively higher redshift would be scruti- 
nized to locate that epoch when the IGM was neutral. 

For an optical spectrum, redshifted to reveal that portion of the rest-frame 
UV shortward of the Lyman a emission line of the QSO itself (A 1216 A ), the 
relative transmission T is defined as 

T = !{\)l icout 

where fcont represents the continuum radiation from the QSOs. The trans- 
mission is reduced by Lyman a absorption in any foreground (lower redshift) 
clouds of neutral hydrogen whose Gunn-Peterson opacity is then 

TCP = -InT 

The first 'complete troughs' in the absorption line spectra of distant QSO 
were presented by Becker et al (2001) and Djorgovski et al (2001) and a more 
comprehensive sampling of 11 SDSS QSOs was presented by Fan et al (2003). 
Recently, an analysis of 19 5.74< z <6.42 QSOs was presented by Fan et al 
(2006b). 

Figure 35 illustrates how absorption structures along the same line of sight 
can be independently probed using Lyman a and higher order lines such as 
Lyman p. 

To understand how this is efi'ective, in a uniform medium tgp is related 
to the abundance of absorbing neutral hydrogen atoms by the following ex- 
pression: 

2 
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TGP = fa Aa H~'^{z) UrI 

nieC 



Observations of the High Redshift Universe 



59 




Fig. 35. The absorption Une spectrum of two SDSS QSOs from the recent study of 

Fan et al (2006b). Here the structures in the Lyman scries absorption lines a and 
P have been aligned in redshift space thereby improving the signal to noise along a 
given sightline. The effects of cosmic variance can clearly be seen by comparing the 
structures seen along the two sightlines. 



where / is the oscillator strength of the Lyman a line, H is the redshift- 
dependent Hubble parameter and Hhi is the neutral number density. 
Numerically, this becomes 

TCP -1.810 /I Qm ^(^)^ — 

where Xhi = (-^-y-^)"^/^ is then the neutral fraction. 

Inspection of this equation is quite revealing. Firstly, even a tiny neutral 
fraction, xhi ^ 10"'*, would give a very deep, seemingly complete GP trough; 
for reference xhi — 1Q~^ today. So clearly the test is not a very sensitive one 
in absolute terms. 

Secondly, since tgp cjc / A, for the same uh, the optical depth in the higher 
order Lyman /3 and 7 lines would be ~6 and 18 times smaller respectively. 

In practice, the above relations are greatly complicated by any dumpiness 
in the medium. This affects our ability to make direct inferences on xhi as 
well as to combine the various Lyman lines into a single test. 
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Instead, workers have examine the relative distribution of tgp with red- 
shift independently from the various Lyman series absorption statistics. An 
increase in tgp with redshift could just be a natural thickening of the Lyman 
absorption forest and, given its weak connection with Xhi, not imply anything 
profound about cosmic rcionization. However, if it can be shown empirically 
that the various diagnostics show a discontinuity in the a;///-redshift trends, 
conceivably we are approaching the neutral era. 

Figure 36 shows that for z <5.5, r oc (-^-g^)"*'^ for both the Lya and 
Ly/3 forests to reasonable precision. However, beyond z ~5.5, both redshift 
trends are much steeper, a (1 + z)^^. The dispersion around the trends also 
increases significantly at higher redshift. Taken together, both results suggest a 
qualitative change in nature of the IGM beyond z ~5.5 (but for an alternative 
explanation see Becker et al 2006). 




Fig. 36. Evolution in the Gunri-Pctcrson optical depth, r, for both the Lyman a 
(left) and /3 (right) forests from the distant SDSS QSO analysis of Fan et al (2006. 
The dotted lines represent fits to the data for Zabs <5.5, beyond which there is 
evidence in both species for an upturn in the opacity of the intergalactic medium. 



Fan et al (2006a,b) discuss several further probes of the nature of the 
IGM at z ~6. One relates to the proximity effect - the region around each 
QSO where it is clear from the spectrum that the IGM is being ionized by 
the QSO itself. Although this region is excluded in the analyses above, the 
extent of this region contains valuable information on the nature of the IGM. 
It appears that the radius of the region affected, R is less at higher redshift 
according to i? oc [(1 + z) Xhi]~^^^ suggesting that the most distant QSOs in 
the sample {z ~6.5) lie in a IGM whose neutral fraction is ~14 times higher 
than those at z ~5.7. 

Another valuable measure is how the regions of complete absorption, the 
so-called 'dark gaps' in the spectrum where transmission is effectively zero, 
are distributed. Fan et al define a 'gap' as a contiguous region in redshift 
space where r >3.5. The distribution of gaps contains some information on 
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the topology of reionization. We would expect regions of high transmission 

to be associated with large HII regions, centered on luminous star-forming 
sources. The dark gaps increase in extent from ~ 10 to 80 comoving Mpc over 
the redshift range samples suggesting the IGM is still not neutral at z ~6.5. 
Although the Gunn-Pcterson and gap statistics make similar statements about 
reionization, suggesting the neutral fraction at z ~6.2 is 1-4%, Fan et al 
consider that beyond z ~6.5, gap statistics will become a more powerful probe. 
This is because the redshift distribution of those few spectroscopic pixels 
where the transmission is non-zero will become the only effective signal. 

5.3 Metallicity of the High Redshift IGM 

A second measure of the redshift range of early star formation is contained 
in the properties of the CIV forest observed in the spectra of high redshift 
QSOs (Songaila 2005, 2006). Carbon is only produced in stellar nuclei (it is 
not produced in the hot Big Bang) and so the ubiquity of CIV along many 
sightlines to z ~5-6 QSOs is a powerful argument for early enrichment. 

CIV was seen in the Lya forest in 1995 (Cowie et al 1995) with N(CIV) /N(HI) 
~10"2 to 10^3 However, it was subsequently seen in even the weakest Lya 
systems (Ellison et al 2000). This is a particularly powerful point since it 
argues that enrichment is not confined to localized regions of high column 
density but is generic to the intergalactic medium as a whole (Figure 37) 

A quantitative interpretation of the CIV abundance, in terms of how much 
early star formation occurred earlier than the highest redshift probed, relies on 
locating a 'floor' in the abundance-redshift relation. Unfortunately, the actual 
observed trend, measured via the contribution of the ion to the mass density 
fl{CIV) from z ~5 to 2, does not seem to behave in the manner expected. For 
example, there is no strong rise in the CIV abundance to lower redshift de- 
spite the obvious continued star formation that occurred within these epochs 
(Songaila 2006, Figure 38). This is a major puzzle (c.f. Oppenheimer et al 
2006). 

5.4 Linear Polarization in the WMAP Data 

In 2003, the WMAP team (Kogut et al 2003) presented the temperature- 
polarization cross angular power spectrum from the first year's data and 
located a 4(7 non-zero signal at very low multipoles (/ <8) which they in- 
terpreted in terms of foreground electron scattering of microwave background 
photons with an optical depth Tg = 0.17 ± 0.04 corresponding to ionized 
structures at -2reion —20 (Figure 39a). The inferred redshift range depends 
sensitively on the history of the reionization process. Bennett et al (2003) ar- 
gued that if reionization occurred instantaneously it corresponds to a redshift 
Zreion — 17 ± 5, whereas adopting a more reasonable Press-Schechter formal- 
ism and an illustrative cooling and enrichment model, Fukugita & Kawasaki 
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Fig. 37. (Left) Keck ESI absorption line spectrum of the z=4.5 QSO BR2237-0607 
from the study of Songaila (2005). A sparsely populated CIV forest from 7500- 
8500A accompanies the dense Lya forest seen below 6800 A . (Right) Distribution 
of column densities of CIV absorbers per unit redshift interval in Q1422+231 from 
the survey of Ellison et al (2000). 



(2003) demonstrated that the same signal can be interpreted with a delayed 
reionization occurring at Zreion —9-10. 

Just before the Saas-Fee lectures, the long-awaited third year WMAP data 
was published (Spergel et al 2006). A refined analysis significantly lowered 
both the normalization of the dark matter power spectrum to fig = 0.74 ± 
0.05, and the optical depth to electron scattering to — 0.09 ± 0.03. The 
same model of instantaneous reionization reduces the corresponding redshift 
to Zreion = 11 ± 3(2(t) - a significant shift from the 1 year data. 

To illustrate the uncertainties, Spergel et al introduce a more realistic 
history of reionization via the ionization fraction Xe- Suppose above Zreion, 
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Fig. 38. Modest evolution in the contribution of intergalactic CIV and SIV over 
2< z <5 as measured in terms of the ionic contribution to the mass density, Q 
(Songaila 2005). 



a;e = and below = \. Then suppose Zj-eion is defined as that in- 

termediate point when = for 7< z < Zreion- Figure 39b illustrates the 
remarkable insensitivity of zreion to the adopted value of Xg for <0.5. De- 
spite improved data, the redshift range implied by the WMAP data spans the 
full range 10< z <20. 

5.5 Stellar Mass Density at z ~5-6 

Neither the Gunn-Peterson test nor the WMAP polarization data necessarily 
demonstrate that reionization was caused by early star-forming sources; both 
only provide constraints on when the intergalactic medium was first reionized. 
The CIV test is a valuable complement since it provides a measure of early 
enrichment which can only come from star-forming sources. Unfortunately, 
powerful though the ubiquitous presence of CIV is in this context, as we 
have seen, quantitative constraints are hard to derive. We thus seek a further 
constraint on the amount of early star formation that might have occurred. 

In Lecture 4 we introduced the techniques astronomers are now using to 
derive stellar masses for distant galaxies. Although the techniques remain 
approximate, it must follow that the stellar mass density at a given epoch 
represents the integral over time and volume of the past star formation. In- 
deed, we already saw a successful application of this in reconciling past star 
formation with the local stellar mass density observed by 2dF (Figure 21). 
Specifically, at a particular redshift, z 



64 Richard S Ellis 




Multipole i ^reac 



Fig. 39. (Left) The 4(t detection of reionization via an excess signal at large scales 
in the angular cross correlation power spectrum of the temperature and polariza- 
tion data in the first year WMAP data (Kogut et al 2003). (Right) Constraints on 
the redshift of reionization, Zreion, from the third year WMAP data (Spergel et al 
2006). The contours illustrate how the Zreion inferred from the lowered optical depth 
depends on the history of the ionized fraction Xe{z), see text for details. 



= / p4z)dV{z) 

J oo 

Using the techniques described in Lecture 4, stellar mass estimates have 
become available for some very high redshift galaxies detected by Spitzer 
(Eyles et al 2005, Mobasher et al 2005, Yan et al 2005). For the most lumi- 
nous Lyman 'dropouts', these estimates are quite substantial, some exceeding 
10^^ Mq implying much earlier activity. Recently, several groups (Stark et al 
2006a, Yan et al 2006, Eyles et al 2006) have been motivated to provide the 
first crude estimates on the volume averaged stellar mass density at these 
early epochs. Part of this motivation is to check whether the massive galaxies 
seen at such high redshift can be reconciled with hierarchical theory, but as 
Stark & EUis (2005) proposed, the established stellar mass can also be used 
to probe earlier star formation and its likely impact upon cosmic reionization. 

A very relevant question is whether the observed mass density at z ~5-6 
is greater than can be accounted for by the observed previous star formation 
history. We will review the rather uncertain data on the star formation density 
p*{z) beyond z ~6 in the next Lecture. However, Stark et al (2006a) find that 
even taking a reasonably optimistic measure of Pi,{z) from recent compilations 
by Bouwens et al (2006) and Bunker et al (2006), it is hard to account for the 
stellar mass density at z ~5 (Figure 40). 

There are currently two major limitations in this comparison. First, most 
of the V- and i-drops are located photometrically; even a small degree of con- 
tamination from lower redshift galaxies could upward bias the stellar mass 
density. On the other hand only star-forming galaxies are located by the Ly- 
man break technique so this bias could easily be offset if there are systems in 
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Fig. 40. A comparison of the assembly history of stellar mass inferred from the 
observed decline in star formation history to z ~10 (solid line) with extant data 
on the stellar mass density at z ~5 and 6 (data points from Stark et al 2006a, 
Yan et al 2006, Eyles et al 2006). Different estimates at a given redshift represent 
lower limits based on spectroscopically-confirmed and photometric redshift samples. 
The red line shows the growth in stellar mass expected from the presently-observed 
luminous star forming galaxies; a shortfall is observed. The blue dotted line shows 
the improvement possible when a dominant component of high z lower luminosity 
systems is included. 



a quiescent state as evidenced by the prominent Balmer breaks seen in many 
of the Spitzer-detected sources (Eyles et al 2005, 2006). This limitation will 
ultimately be overcome with more careful selection methods and deeper spec- 
troscopy. Secondly, and more profoundly, the precision of the stellar masses 
may not be up to this comparison. Much has to be assumed about the nature 
of the stellar populations involved which may, quite reasonably, be somewhat 
different from those studied locally. The discrepancy noted by Stark et al is 
only a factor of x2-3, possibly within the range of uncertainty. 

Regardless, if this mismatch is reinforced by better data, the implications 
are very interesting in the context of reionization. It could mean early star- 
forming systems are extincted, lie beyond z ~10 where current searches end, or 
perhaps most likely that early star formation is dominated by lower luminosity 
systems (Figure 40). By refining this technique and using diagnostics such as 
the strength of the tell-tale Balmer break, it may ultimately be possible to 
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age-date the earlier activity and compare its efficacy with that required to 
reionize the Universe. 

5.6 Lecture Summary 

In this lecture we have introduced four very different and independent probes 
of cosmic reionization, each of which suggests star formation activity may 
extend well into the redshift range 6 < z < 20. Two of these probes rely on 
a contribution from early star formation (the mctallicity of the intergalactic 
medium and the assembled stellar mass density at z ~5-6). 

The earliest result was the presence of neutral hydrogen troughs in the 
spectra of distant QSOs. Although the arguments for reionization ending at 
z seem compelling at first sight, they ultimately rely on an empirically- 
deduced transition in the changes in the opacity of the Lya and Ly/3 line 
below and above z ~5.5. 

The second result - the ubiquity of carbon in even the weakest absorbing 
clouds at z ~5 is firm evidence for early star formation. However, it seems hard 
to locate the high redshift 'abundance floor' and hence to quantify whether 
this early activity is sufficient for reionization. Indeed, a major puzzle is the 
lack of growth in the carbon abundance over the redshift range where galaxies 
are assembling the bulk of their stars. 

The WMAP polarizations results have received the most attention, mainly 
because the first year data indic;atcd a surprisingly large optical depth and a 
high redshift for reionization. However, there were some technical limitations 
in the original analysis and it now seems clear that the constraints on the 
redshift range when the foreground polarization is produced are not very 
tight. 

Finally, an emerging and very promising technique is simply the census of 
early star formation activity as probed by the stellar masses (and ages) of the 
most luminous dropouts at z ~5-6. Although significant uncertaintes remain, 
the prospects for improving these constraints are good and, at this moment, 
it seems there must have been quite a significant amount of early (z >6) star 
formation activity, quite possibly in low luminosity precursors. 
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6 Into the Dark Ages: Lyman Dropouts 
6.1 Motivation 

Surveys of galaxies at and beyond a redshift z ~6 represent the current ob- 
servational frontier. We are motivated to search to conduct a census of the 
earliest galaxies seen 1 Gyr after the Big Bang as well as to evaluate the con- 
tribution of early star formation to cosmic reionization. Although impressive 
future facilities such as the next generation of extremely large telescop^l and 
the James Webb Space Telescope (Gardner et al 2006) are destined to address 
these issues in considerable detail, any information we can glean on the abun- 
dance, luminosity and characteristics of distant sources will assist in planning 
their effective use. 

In this lecture and the next, we will review the current optical and near- 
infrared techniques for surveying this largely uncharted region. They include 

• Lyman dropouts: photometric searches based on locating the rest- frame 
ultraviolet continuum of star-forming sources introduced in Lecture 3. The 
key issue here is reducing contamination from foreground sources since 
most sources selected via this technique are too faint for confirmatory 
spectroscopy. 

• Lyman alpha emitters: spectroscopic or narrow-band searches for sources 
with intense Lya emission. As the line is resonantly-scattered by neutral 
hydrogen its profile and strength gives additional information on the state 
of the high redshift intergalactic medium. 

• Strong gravitational lensing: by coupling both above techniques with the 
magnification afforded by lensing clusters, it is possible to search for lower 
luminosity sources at high redshift. Since the magnified areas are small, 
the technique is only advantageous if the luminosity function has a steep 
faint end slope. 

For the Lyman dropouts discussed here, as introduced in the last lecture, 
there is an increasingly important role played by the Spitzer Space Telescope 
in estimating stellar masses and earlier star formation histories. 

The key questions we will address in this lecture focus on the (somewhat 
controversial) conclusions drawn from the analyses of Lyman dropouts thus 
far, namely: 

1. How effective are the various high z selection methods? The characteris- 
tic luminosity L* at z ~6 corresponds to iab —26 where spectroscopic 
samples are inevitably biased to those with prominent Lya emission. Ac- 
cordingly there is great reliance on photometric redshifts and a real danger 
of substantial contamination by foreground red galaxies and Galactic cool 
stars. 



e.g. The US-Canadian Thirty Meter Telescope - http//www. tmt.org 
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2. Is there a decline in the UV luminosity density, puv-, over the range 3< 
z <6? The early results were in some disagreement. Key issues relate to 
the degree of foreground contamination and cosmic variance in the very 
small deep fields being examined. 

3. Is the observed UV density, pfjy, at z ci6 sufficient to account for reion- 
ization? The answer depends on the contribution from the faint end of the 
luminosity function and whether the UV continuum slope is steeper than 
predicted for a normal solar-metallicity population. 

4. Significant stellar masses have been determined for several z ~6 galaxies. 
Are these in conflict with hierarchical structure formation models? 

6.2 Contamination in z ~6 Dropout Samples 

The traditional dropout technique exploited very effectively at z ~3 (Lecture 
3) is poorly suited for z ~6 samples because the use of a simple i—z > 1.5 color 
cut still permits significant contamination by passive galaxies at z ~2 and 
Galactic stars. The addition of an optical-infrared color allows some measure 
of discrimination (Stanway et al 2005) since a passive z ~2 galaxy will be 
red over a wide range in wavelength, whereas a star-forming z ~6 galaxy 
should be relatively blue in the optical-infrared color corresponding to its rest- 
frame ultraviolet (Figure 41). Application of this two color technique suggests 
contamination by foreground galaxies is ~10% at the bright end (zab <25.6) 
but negligible at the UDF limit {zab <28.5) 

Unfortunately, the spectral properties of cool Galactic L dwarfs are dom- 
inated by prominent molecular bands rather than simply by their effective 
temperature. This means that they cannot be separated from z ~6 galaxies 
in a similar color-color diagram. Indeed, annoyingly, these dwarfs occupy pre- 
cisely the location of the wanted z ~6 galaxies (Figure 42)! The only practical 
way to discriminate L dwarfs is either via spectroscopy or their unresolved 
nature in ACS images. 

Stanway et al (2004) conducted the first comprehensive spectroscopic and 
ACS imaging survey of a GOODS i-drop sample limited at zab <25.6, finding 
that stellar contamination at the bright end of the luminosity function of a 
traditional (zab — zab >1.5) color cut could be as high as 30-40%. Unfortu- 
nately, even with substantial 6-8 hour integrations on the Keck telescope, red- 
shift verification of the distant population was only possible in those dropout 
candidates with Lya emission. Stanway et al (2005) subsequently analyzed the 
ACS imaging properties of a fainter subset arguing that stellar contamination 
decreases with increasing apparent magnitude. 

Further progress has been possible via the use of the ACS grism on board 
HST (Malhotra et al 2005). As the OH background is eliminated in space, 
despite its low resolution, it is possible even in the fairly low signal/noise data 
achievable with the modest 2.5m aperture of HST to separate a Lyman break 
from a stellar molecular band. It is claimed that of 29 Zy\B <27.5 candidates 
with (i-z)>0.9, only 6 are likely to be low redshift interlopers. 
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Fig. 41. The combination of a i — z and z — J color cut permits the distinction of 
z ~5.7-6.5 star forming and z ~2 passive galaxies. Both may satisfy the i — z >1.5 
dropout selector, but the former should lie blueward of the z — ,7=1.0 divider, whereas 
z ~2 are red in both colors. Crosses represent the location of candidates in the 
GOODS field and model tracks illlustrate the predicted colors for typical SEDs 
observed at the respective redshifts (Stanway et al 2005). 

Regrettably, as a result of these difHculties, it has become routine to rely 
entirely on photometric and angular size information without questioning fur- 
ther the degree of contamination. This is likely one reason why there remain 
significant discrepancies between independent assessments of the abundance 
of z ~6 galaxies (Giavalisco et al 2004, Bouwens et al 2004, Bunker et al 2004). 
Although there are indications from the tests of Stanway et al (2004, 2005) 
and Malhotra et al (2005) that contamination is significant only at the bright 
end, the lack of a comprehensive understanding of stellar and foreground con- 
tamination remains a major uncertainty. 

6.3 Cosmic Variance 

The deepest data that has been searched for i-band dropouts includes the two 
GOODS fields (Dickinson et al 2004) and the Hubble Ultra Deep Field (UDF, 
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Fig. 42. (Left) Keck spectroscopic verification of two contaminating L dwarfs lying 
within the GOODS i — z dropout sample but pinpointed as likely to be stellar from 
ACS imaging data. The smoothed spectra represent high signal to noise brighter 
examples for comparison purposes. Strong molecular bands clearly mimic the Ly- 
man dropout signature. (Right) Optical-infrared color diagram with the dropout 
color selector, zas — zab > 1-5, shown as the vertical dotted line. Bright L dwarfs 
(lozenges) frustratingly occupy a similar region of color space as the z ~6 candidates 
(points with error bars) (Stanway et al 2004). 

Beckwith ct al 2006). As these represent publicly-available fields they have 
been analyzed by many groups to various flux limits. The Bunker/Stanway 
team probed the GOODS fields to z^b=25.6 (spectroscopically) and 27.0 
(photometrically, and the UDF to .2ab=28.5. At these limits, it is instructive 
to consider the comoving cosmic volumes available in each field within the 
redshift range selected by the typical dropout criteria. For both GOODS-N/S 
fields, the total volume is ~5. 10^ Mpc^, whereas for the UDF it is only 2.6 
10^ Mpc^. These contrast with 10^ Mpc^ for a single deep pointing taken with 
the SuPrime Camera on the Subaru 8m telescope. 

Somerville et al (2004) present a formalism for estimating, for any popula- 
tion, the fractional uncertainty in the inferred number density from a survey 
of finite volume and angular extent. When the clustering signal is measurable, 
the cosmic variance can readily be calculated analytically. However, for fron- 
tier studies such as the i-dropouts, this is not the case. Here Somerville et al 
propose to estimate cosmic variance by appealing to the likely halo abundance 
for the given observed density using this to predict the clustering according 
to CDM models. In this way, the uncertainties in the inferred abundance of 
i-dropouts in the combined GOODS fields could be ~20-25% whereas that in 
the UDF could be as high as 40-50%. 
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It seems these estimates of cosmic variance can only be strict lower limits 
to the actual fluctuations since Somerville et al make the assumption that 
halos containing star forming sources are visible at all times. If, for example, 
there is intermittent activity with some duty cycle whose "on/ofF' fraction is 
/, the cosmic variance will be underestimated by that factor / (Stark et al, 
in prep). 

6.4 Evolution in the UV Luminosity Density 3< z <10? 

The complementary survey depths means that combined studies of GOODS 
and UDF have been very effective in probing the shape of the UV luminosity 
function (LF) at z ^qj- Even so, there has been a surprising variation in 
the derived faint end slope a. Bunker et al (2004) claim their data (54 i- 
dropouts) is consistent with the modestly-steep a=-1.6 found in the z ~3 
Lyman break samples (Steidel et 2003), whereas Yan & Windhorst (2004) 
extend the UDF counts to 0yiB=3O.O and, based on 108 candidates, find a=- 
1.9, a value close to a divergent function! Issues of sample completeness are 
central to understanding whether the LF is this steep. 

In a comprehensive analysis based on all the extant deep data, Bouwens et 
al (2006) have attempted to summarize the decline in rest-frame UV luminos- 
ity density over 3< z <10 as a function of luminosity (Figure 43). They at- 
tribute the earlier discrepancies noted between Giavalisco et al (2004) , Bunker 
et al (2004) and Bouwens et al (2004) to a mixture of cosmic variance and dif- 
ferences in contamination and photometric selection. Interestingly, they claim 
a luminosity-dependent trend in the sense that the bulk of the decline occurs 
in the abundance of luminous dropouts, which they attribute to hierarchical 
growth. 

A similar trend is seen in ground-based data obtained with Subaru. Al- 
though HST offers superior photometry and resolution which is effective in 
eliminating stellar contamination, the prime focus imager on Subaru has a 
much larger field of view so that each deep exposure covers a field twice as 
large as both GOODS N-t-S. As they do not have access to ACS data over such 
wide fields, the Japanese astronomers have approached the question of stellar 
contamination in an imaginative way. Shioya et al (2005) used two interme- 
diate band filters at 709nm and 826nm to estimate stellar contamination in 
both z ~5 and z ~6 broad-band dropout samples. By considering the slope 
of the continuum inbetween these two intermediate bands, in addition to a 
standard i — z criterion, they claim an ability to separate L and T dwards. 
In a similar, but independent, study, Shimasaku et al (2005) split the z band 
into two intermediate filters thereby measuring the rest-frame UV slope just 
redward of the Lyman discontinuity. These studies confirm both the redshift 
decline and, to a lesser extent, the luminosity-dependent trends seen in the 
HST data. 

In this section we will only refer to the observed (extincted) LPs and luminosity 
density 
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Fig. 43. Evolution in the rest-frame UV (1350 A ) luminosity density (right ordi- 
nate) and inferred star formation rate density ignoring extinction (left ordinate) for 
drop-out samples in two luminosity ranges from the compilation by Bouwens et al 
(2006). A marked decline is seen over 3 < z < 6 in the contribution of luminous 
sources. 



Although it seems there is a 5 x abundance decline in luminous UV emit- 
ting galaxies from z ~3 to 6, it's worth noting again that the relevant counts 
refer to sources uncorrected for extinction. This is appropriate in evaluating 
the contribution of UV sources to the reionization process but not equivalent, 
necessarily, to a decline in the star formation rate density. Moreover, although 
the luminosity dependence seems similar in both ground and HST-based sam- 
ples, it remains controversial (e.g. Beckvifith et al 2006). 

6.5 The Abundance of Star Forming Sources Necessary for 
Reionization 

Have enough UV-emitting sources been found at z ~6-10 to account for cos- 
mic reionization? Notwithstanding the observational uncertainties evident in 
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Figure 43, this has not prevented many teams from addressing this important 
question. The main difficulty lies in understanding the physical properties of 
the sources in question. The plain fact is that we cannot predict, sufficiently 
accurately, the UV luminosity density that is sufficient for reionization! 

Some years ago, Madau, Haardt & Rccs (1999) estimated the star for- 
mation rate density based on simple parameterized assumptions concern- 
ing the stellar IMF and/or metallicity Z essential for converting a 1350 
A luminosity into the integrated UV output, the fraction /esc of escap- 
ing UV photons, the dumpiness of the surrounding intergalactic hydrogen, 
C =< p\ii > I < Phi and the temperature of the intergalactic medium 
TiGM- In general terms, for reionization peaking at a redshift Zreion-, the nec- 
essary density of sources goes as: 

p^f-lC{l + zf{QBh''fMpc-^ 

For likely ranges in each of these parameters, StiavcUi ct al (2004) tabulate 
the required source surface density which, generally speaking, lie above those 
observed at z ~6 (e.g. Bunker et al 2004). 

It is certainly possible to reconcile the end of the reionization at z ~6 with 
this low density of sources (Figure 43) by appealing to cosmic variance, a low 
metallicity and/or top-heavy IMF (Stiavelli et al 2005) or a steep faint end 
slope of the luminosity function (Yan & Windhorst 2004) but none of these 
arguments is convincing without further proof. As we will see, the most logical 
way to proceed is to explore both the extent of earlier star formation from 
the mass assembled at z ~5-6 (Lecture 5) and to directly measure, if possible, 
the abundance of low luminosity sources at higher redshift. 

6.6 The Spitzer Space Telescope Revolution: Stellar Masses at 
z ~6 

One of the most remarkable aspects of our search for the most distant and 
early landmarks in cosmic history is that a modest cooled 85cm telescope, the 
Spitzer Space Telescope, can not only assist but provide crucial diagnostic 
data! The key instrument is the InfraRed Array Camera (IRAC) which offers 
four channels at 3.6, 4.5, 5.8 and 8/im corresponding to the rest-frame optical 
0.5-1/um at redshifts z ~6-7. In the space of only a year, the subject has 
progressed from the determination of stellar masses for a few z ^5-7 sources 
to mass densities and direct constraints on the amount of early activity, as 
discussed in Lecture 5. 

An early demonstration of the promise of IRAC in this area was provided 
by Eyles et al (2005) who detected two spectroscopically-confirmed z ~5.8 
«-band dropouts at 3.6/Ltm, demonstrating the presence of a strong Balmer 
break in their spectral energy distributions (Figure 44). In these sources, the 
optical detection of Lya emission provides an estimate of the current ongoing 
star formation rate, whereas the flux longward of the Balmer break provides 
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a measure of the past averaged activity. The combination gives a measure 
of the luminosity weAghted age of the stellar population. In general terms, 
a Balmer break appears in stars whose age cannot, even in short burst of 
activity, be younger than 100 Myr. Eyles et al showed such systems could 
well be much older (250-650 Myr) depending on the assumed form of the past 
activity. As the Universe is only 1 Gyr old at z ~6, the IRAC detections gave 
the first indirect glimpse of significant earlier star formation - a glimpse that 
was elusive with direct searches at the time. 




Fig. 44. (Left) Detection of a spectroscopically-confirmed i-drop at z=5.83 from 
the analysis of Eyles et al (2005). (Right) Spectral energy distribution of the same 
source. Data points refer to IRAC at 3.6 and 4.5/xm, VLT (K) and HST NICMOS 
(J,H) overplotted on a synthesised spectrum; note the prominent Balmer break. 
Synthesis models indicate the Balmer break takes 100 Myr to establish. However, 
the luminosity-weighted age could be significantly older depending on the assumed 
past star formation rate. In the example shown, a dominant 450 Myr component 
{zf ~10) is rejuvenated with a more recent secondary burst whose ongoing star 
formation rate is consistent with the Lya flux observed in the source. 

Independent confirmation of both the high stellar masses and prominent 
Balmer breaks was provided by the analysis of Yan et al (2005) who studied 
3 z ~5.9 sources. Moreover, Yan et al also showed several objects had {z — 
J) colors bluer than the predictions of the Bruzual-Charlot models for all 
reasonable model choices - a point first noted by Stanway et al (2004). 

Eyles et al and Yan et al proposed the presence of established stellar 
populations in z ~6 i-drops and also to highlight the high stellar masses 
(M — A IO'^^Mq) they derived. At first sight, the presence of z ~6 sources 
as massive as the Milky Way seems a surprising result. Yan et al discuss the 
question in some detail and conclude the abundance of such massive objects is 
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not inconsistent with hierarchical theory. In actuaUty it is hard to be sure be- 
cause cosmic variance permits a huge range in the derived vohime density and 
theory predicts the halo abundance (e.g. Barkana & Loeb 2000) rather than 
the stellar mass density. To convert one into the other requires a knowledge 
of the star formation efficiency and its associated duty-cycle. 

One early UDF source detected by IRAC has been a particular source of 
puzzlement. Mobasher et al (2005) found a J-dropout candidate with a promi- 
nent detection in all 4 IRAC bandpasses. Its photometric redshift was claimed 
to be z ~6.5 on the basis of both a Balmer and a Lyman break. However, de- 
spite exhaustive efforts, its redshift has not been confirmed spectroscopically. 
The inferred stellar mass is 2-7 VS^'^Mq, almost an order of magnitude larger 
than the spectroscopically-confirmed sources studied by Eyles et al and Yan 
et al. If this source is truly at z ~6.5, finding such a massive galaxy whose 
star formation likely peaked before z ~9 is very surprising in the context of 
contemporary hierarchical models. Such sources should be extremely rare so 
finding one in the tiny area of the UDF is all the more puzzling. Dunlop et 
al (2006) have proposed the source must be foregroimd both on account of 
an ambiguity in the photometric redshift determination and the absence of 
similarly massive sources in a panoramic survey being conducted at UKIRT 
(McClurc et al 2006). 

This year, the first estimates of the stellar mass density at z ~5-6 have 
been derived (Yan et al 2006, Stark et al 2006a, Eyles et al 2006). Although 
the independenty-derived results are consistent, both with one another and 
with lower redshift estimates (Figure 45) the uncertainties are considerable as 
discussed briefly in the previous lecture. There are four major challenges to 
undertaking a census of the star formation at early times. 

Foremost, the bulk of the faint sources only have photometric redshifts. 
Even a small amount of contamination from foreground sources would skew 
the derived stellar mass density upward. Increasing the spectroscopic coverage 
would be a big step forward in improving the estimates. 

Secondly, IRAC suffers from image confusion given its lower angular res- 
olution than HST (4 arcsec c.f. 0.1 arcsec). Accordingly, the IRAC fluxes 
cannot be reliably estimates for blended sources. Stark et al address this by 
measuring the masses only for those uncontaminated, isolated sources, scaling 
up their total by the fraction omitted. This assumes confused sources are no 
more or less likely to be a high redshift. 

Thirdly, as only star forming sources are selected using the v— and i— 
dropout technique, if star formation is episodic, it is very likely that quiescent 
sources are present and thus the present mass densities represent lower limits. 
The missing fraction is anyone's guess. As we saw at z ~2, the factor could 
be as high as x2. 

Finally, as with all stellar mass determinations, many assumptions are 
made about the nature of the stellar populations involved and their star for- 
mation histories. Until individual z ~5-6 sources can be studied in more detail, 
perhaps via the location of one or two strongly lensed examples, or via future 



76 Richard S Ellis 



10'^ 



IT 



I 

Iw. 

iS 

4> - 



10^ 



10*' r 




irMudnick et al. (^003) 
V Stark et al. (2006) 
oYan etal. (2006) 
^ Cole et al . (2001 ) ^Thfs work 

.jr.4 i l,V-t„. j . I:>-.:j,-,> .. , ■ I I , 



1 2 3 4 5 6 
Redshift 

Fig. 45. Evolution in the comoving stellar mass density from the compilation de- 
rived by Eyles et al (2006). The recent z ~5-6 estiimates constitute lower limits 
given the likelihood of quiescent sources missed by the drop-out selection technique. 
Results at z ~6 are offset slightly in redshift for clarity. 



more powerful facilities, this will regrettably remain the situation. At present, 
such density estimates are unlikely to be accurate to better than a factor of 
2. Even so, they provide good evidence for significant earlier star formation 
(Stark et al 2006a, Lecture 5). 

6.7 Lecture Summary 

In this lecture we have discussed the great progress made in using v,i,z and 
J band Lyman dropouts to probe the abundance of star forming galaxies over 
3< z <10. At redshifts z ~6 alone, Bouwens et al (2006) discuss the properties 
of a catalog of 506 sources to zab=^Q-5. 

In practice, the good statistics are tempered by uncertain contamination 
from foreground cool stars and dusty or passively-evolving red z ~2 galaxies 
and the vagaries of cosmic variance in the small fields studied. It may be 
that we will not overcome these difficulties until we have larger ground-based 
telescopes. 

Nonetheless, from the evidence at hand, it seems that the comoving UV 
luminosity density declines from z ~3 to 10, and that only by appealing to 
special circumstances can the low abundance of star forming galax;ies at z >6 
be reconciled with that necessary to reionize the Universe. 
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One obvious caveat is our poor knowledge of the contribution from lower 
luminosity systems. Some authors (Yan & Windhorst 2004, Bouwens et al 
2006) have suggested a steepening of the luminosity function at higher red- 
shift. Testing this assumption with lensed searches is the subject of our next 
lecture. 

Finally, we have seen the successful emergence of the Spitzer Space Tele- 
scope as an important tool in confirming the need for star formation at z >6. 
Large numbers of z ~5-6 galaxies have now been detected by IRAC. The 
prominent Balmer breaks and high stellar masses argue for much earlier ac- 
tivity. Reconciling the present of mature galaxies at z ~6 with the absence of 
significant star formation beyond, is one of the most interesting challenges at 
the present time. 
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7 Lyman Alpha Emitters and Gravitational Lensing 

7.1 Strong Gravitational Lensing - A Primer 

Slowly during the twentieth century, gravitational lensing moved from a cu- 
riosity associated with the verification of General Relativity (Eddington f 9f 9) 
to a practical tool of cosmologists and those studying distant galaxies. There 
are many excellent reviews of both the pedagogical aspects of lensing (Bland- 
ford & Narayan 1998, MeUier 2000, Refregier 2002) and a previous Saas-Fee 
contributor (Schneider 2006). 

To explore the distant Universe, we are primarily concerned with strong 
lensing - where the lens has a projected mass density, above a critical value, 
Scrit, so that multiple images and high source magnifications are possible. 
For a simple thin lens 



where D represents the angular diameter distance and the subscripts 
0,L,S refer to the observer, source and lens, respectively. Rather conve- 
niently, for a lens at z ~0.5 and a source at z >2, the critical projected 
density is about 1 g cm~^ - a value readily exceeded by most massive clusters. 
The merits of exploring the distant Universe by imaging through clusters was 
sketched in a remarkably prophetic article by Zwicky (1937). 

In lensing theory it is convenient to introduce a source plane, the true 
sky, and an image plane, the detector at our telescope, where the multiple 
images are seen. The relationship between the two is then a mapping trans- 
formation which depends on the relative distances (above). Crucially, what 
the observer sees depends on the degree of alignment between the source and 
lens as illustrated in Figure 46. 

An elliptical lens with S > Scrit for a given source and lens distance 
produces a pair of critical lines in the image plane where the multiple images 
lie. These lines map to caustics in the source plane. The outer critical line is 
equivalent to the Einstein radius 9e 



and, for a given source and lens, is governed by the enclosed mass M. The 
location of the inner critical line depends on the gradient of the gravitational 
potential (Sand et al 2005). 

The critical lines are important because they represent areas of sky where 
very high magnifications can be encountered - as high as x30! For ~20 well- 
studied clusters the location of these lines can be precisely determined for a 
given source redshift. Accordingly, it is practical to survey just those areas to 
secure a glimpse of otherwise inaccessibly faint sources boosted into view. The 
drawback is that, as in an optical lens, the sky area is similarly magnified, so 
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the surface density of faint sources must be very large to yield any results. 
Regions where the magnification exceed x 10 are typically only 0.1-0.3 arcmin^ 
per cluster in extent in the image plane and inconveniently shaped for most 
instruments (Figure 47). The sampled area in the source plane is then ten 
times smaller so to see even one magnified source/cluster requires a surface 
density of distant sources of ^50 arcmin^^. 

Two other applications are particularly useful in faint galaxy studies. 
Firstly, strongly magnified systems at z ~2-3 can provide remarkable insight 
into an already studied population by providing an apparently bright galaxy 
which is brought within reach of superior instrumentation. cB58, a Lyman 
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break galaxy at z=2.72 boosted by x30 to V=20.6 (Yee et al 1996, Seitz et al 
1998) was the first distant galaxy to be studied with an echellette spectrograph 
(Pettini et al 2002), yielding chemical abundances and outflow dynamics of 
unprecedented precision. 

More generally, a cluster can magnify a larger area of ~2-4 arcmin^ by a 
modest factor, say x3-5. This has been effective in probing sub- mm source 
counts to the faintest possible limits (Small et al 1997) and the method shows 
promise for similar extensions with the IRAC camera onboard Spitzer. 
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Fig. 47. Hubble Space Telescope image of the rich cluster Abell 1689 with the 
critical lines for a source at z ~5 overlaid in yellow. The narrow regions inbetween 
the pairs of yellow lines refer to regions where the magnification exceeds x 10. 



7.2 Creating a Cluster Mass Model 

In the applications discussed above, in order to analyze the results, the inferred 
magnification clearly has to be determined. This will vary as a function of 
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position in the cluster image and the relative distances of source and lens. The 
magnification follows from the construction of a mass model for the cluster. 

The precepts for this method are discussed in the detailed analysis of the 
remarkable image of Abell 2218 taken with the WFPC-2 camera onboard HST 
in 1995 (Kneib et al 1996). An earlier image of AC 114 showed the important 
role HST would play in the recognition of multiple images (Smail et al 1995). 
Prior to HST, multiple images could only be located by searching for systems 
with similar colors, using the fact that lensing is an achromatic phenomenon. 
HST revealed that morphology is a valuable additional identifier; the improved 
resolution also reveals the local shear (see Figure 48). 




Fig. 48. Hubble Space Telescope study of the rich cluster AC114 (Smail et al 
2005) (Top) Morphological recognition of a triply-imaged source. The lower inset 
panels zoom in on each of 3 images of the same source. (Bottom) Construction of 
mass contours (red lines) and associated shear (red vectors) from the geometrical 
arrangement of further multiple images labelled Al-3, Bl-4, Cl-3, Ql-3. 



Today, various approaches are possible for constructing precise mass mod- 
els for lensing clusters (Kneib et al 1996, Julio et al 2006, Broadhurst et al 
2005). These are generally based on utilizing the geometrical positions of sets 
of multiply-imaged systems whose redshift is known or assumed. This then 
maps the form and diameter of the critical line for a given z. Spectroscopic 
redshifts are particularly advantageous, as are pairs that straddle the criti- 
cal line whose location can then be very precisely pinpointed. A particular 
mass model can be validated by 'inverting' the technique and predicting the 
redshifts of other pairs prior to subsequent spectroscopy (Ebbels et al 1999). 

The main debate among cognescenti in this area lies in the extent to 
which one should adopt a parametric approach to fitting the mass distribution, 
particularly in relation to the incorporation of mass clumps associated with 
individual cluster galaxy halos (Broadhurst et al 2005). Stark et al (2006b) 
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discuss the likely uncertainties in the mass modeling process arising from the 
various techniques. 

7.3 Lensing in Action: Some High z Examples 

Before turning to Lyman a emitters (lensed and unlensed), we will briefly 
discuss what has been learned from strongly-lensed dropouts. 

Figure 49 shows a lensed pair in the cluster Abell 2218 (z=0.18) as detected 
by NICMOS onboard HST and the two shortest wavelength channels of IRAC 
(Kneib et al 2004, Egami et al 2005). Although no spectroscopic redshift is 
yet available for this source, three images have been located by HST and their 
arrangement around the well-constrained 2;=6 critical line suggests a source 
beyond z ~6 (Kneib et al 2004). 

As with the unlensed i-band drop out studies by Eyles et al (2005) and 
Yan et al (2005), the prominent IRAC detections (Egami et al 2005) permit an 
improved photometric redshift and important constraints on the stellar mass 
and age. A redshift of z=6.8±0.1 is derived, independently of the geometric 
constraints used by Kneib et al. The stellar mass is ~5-10 10^ Af© and the 
current star formation rate is ~2.6 Mq yr~^. The luminosity-weighted age 
corresponds to anything from 40-450 Myr for a normal IMF depending on 
the star formation history. Interestingly, the derived age for such a prominent 
Balmer break generally exceeds the e-folding timescale of the star formation 
history (Fig. 49) indicating the source would have been more luminous at 
redshifts 1 < z <V1 (unless obscured). 

Given the small search area used to locate this object, such low mass 
sources may be very common. Accordingly, several groups are now survey- 
ing more lensing clusters for further examples of z band dropouts and even 
J-band dropouts (corresponding to z ~8-10) . Richard et al are surveying 6 
clusters with NICMOS and IRAC with deep ground-based K band imaging 
from Subaru and Keck. In these situations, one has to distinguish between 
magnifications of x 5 or so expected across the 2-3 arcmin fields of NICMOS 
and IRAC, and the much larger magnifications possible close to the criti- 
cal lines. Contamination from foreground sources should be similar to what 
is seen in the GOODS surveys discussed in Lecture 6. The discovery of im- 
age pairs in the highly-magnified regions would be a significant step forward 
since spectroscopic confirmation of any sources at the limits being probed 
[Eab =^26.5-27.0) will be exceedingly difficult. 

7.4 Lyman alpha Surveys 

The origin and characteristic properties of the Lyman a emission line has been 
discussed by my colleagues in their lectures (see also Miralda-Escude 1998, 
Haiman 2002, Barkana & Loeb 2004 and Santos 2004). The n=2 to 1 transition 
corresponding to an energy diiference of 10.2 eV and rest- wave-length of A1216 
A typically arises from ionizing photons absorbed by nearby hydrogen gas. The 
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Fig. 49. (Top) Lensed pair of a z—Q.8 source as seen by NICMOS and IRAC in 
the rich cluster Abell 2218 (Egami et al 2005). The pair straddles the critical line 
at z ~6 and a fainter third image at a location predicted by the lensing model has 
been successfully recovered in the HST data (Kneib et al 2004). (Bottom) Spectral 
energy distribution of the source revealing a significant Balmer break and improved 
estimates of the star formation rate, stellar mass and luminosity weighted age. 



line has a number of interesting features which make it particularly well-suited 
for locating early star forming galaxies as well as for characterizing the nature 
of the IGM. 

In searching for distant galaxies, emission lines offer far more contrast 
against the background sky than the faint stellar continuum of a drop-out. A 
line gives a convincing spectroscopic redshift (assuming it is correctly iden- 
tified) and models suggest that as much as 7% of the bolometric output of 
young star-forming region might emerge in this line. For a normal IMF and 
no dust, a source with a star formation rate (SFR) of 1 Mq yr~^ yields an 
emission line luminosity of 1.5 10^^ ergs sec~^. 
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Narrow band imaging techniques (see below) can reach fliuces of < 10"^'' 
cgs in comoving survey volunics of ~10''^ Afpc"^, corresponding to a SFR ~3 
yi'"^ at z ~6. Spectroscopic techniques can probe fainter due to the 
improved contrast. This is particularly so along the critical lines where the 
additional boost of gravitational Icnsing enables fluxes as faint as 3. lO"'^^ 
cgs to be reached (corresponding to SFR ~0.1 Mq yr~^). However, in this 
case the survey volumes are much smaller (~50 Mpc^). In this sense, the two 
t(;c,liiii([ucs (discussed below) arc usefully complementary. 

Having a large dynamic range in surveys for Lycc emission is important not 
just to probe the luminosity function of star-forming galaxies but also because 
it can be used to characterize the IGM. As a resonant transition, foreground 
hydrogen gas clouds can scatter away Lya photons in both direction and 
frequency. In a partially ionized IGM, scattering is maximum at A1216 A in the 
rest-frame of the foreground cloud, thus affecting the blue wing of the observed 
line. However, in a fully neutral IGM, scattering far from resonance can occur 
leading to damping over the entire observed line. Figure 50 illustrates how, in 
a hypothetical situation where the IGM becomes substantially neutral during 
6< z <7, surveys reaching the narrow-band flux limit would still find emitters 
at z=7. Their intense emission would only be partially damped by even a 
neutral medium. However, lines with fluxes at the spectroscopic lensing limit 
would not survive. Accordingly, one possible signature of reionization would 
be a significant change in the shape of the Lya luminosity function at the 
faint end (Furlanetto et al 2005). 

7.5 Results from Narrow Band Lya Surveys 

The most impressive results to date have come from various narrow-band fil- 
ters placed within the SuPrime camera at the prime focus of the Subaru 8m 
telescope (Kodaira et al 2003, Hu et al 2004, Ouchi et al 2005, Taniguchi et 
al 2005, Shimasaku et al 2006, Kashikawa et al 2006, lye et al 2006). Impor- 
tant conclusions have also been deduced from an independent 4m campaign 
(Malhotra & Rhoads 2004). 

Narrow band filters are typically manufactured at wavelengths where the 
night sky spectrum is quiescent, thus maximizing the contrast. These locations 
correspond to redshifts of 2=4.7, 5.7, 6.6 and 6.9 (Fig. 52). A recent triumph 
was the successful recovery of two candidates at z ~6.96 by lye et al (2006). 
Candidates are selected by comparing their narrow band fluxes with that in a 
broader band encompassing the narrow band wavelength range. The contrast 
can then be used as an indicator of line emission (Fig. 53). Spectroscopic 
follow-up is still desirable as the line could arise in a foregound galaxy with 
[O II] 3727 A or [O III] 5007 A emission. The former line is a doublet and 
the latter is part of a pair with a fixed line ratio, separated in the rest-frame 
by only 60 A or so, so these contaminants are readily identified. Furthermore, 
Lya is often revealed by its asymmetric profile (c.f. Fig. 51). 
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Fig. 50. The Lyman a damping wing is absorbed by neutral hydrogen and thus 
can act as a valuable tracer of the nature of the IGM. The simulation demonstrates 
the effect of Ifl damping on emission lines in high mass and low mass systems 
(characteristic of sources detected in typical narrow band and lensed spectroscopic 
surveys respectively) assuming reionization ends inbetween z—6 and 7. The dramatic 
change in visibility of the weaker systems suggests their study with redshift may offer 
a sensitive probe of reionization. Courtesy: Mike Santos. 



Spectroscopic follow-up is obviously time-intensive for a large sample of 
candidates, hundreds of which can now be found with panoramic imagers such 
as SuPrimeCam. Therefore it is worth investigating additional ways of elim- 
inating foreground sources. Tanuguchi et al (2005) combine the narrow band 
criteria adopted in Fig. 53 with a broad-band i — z drop-out signature. Spec- 
troscopic follow-up of candidates located via this double color cut revealed a 
50-70% success rate for locating high z emitters. The drawback is that the 
sources so found cannot easily be compared in number with other, more tradi- 
tional, methods. Nagao et al (2005) used a narrow - broad band color criterion 
in the opposite sense, locating sources with a narrow band depression (rather 
than excess) . Such rare sources are confirmed to be sources with extremely in- 
tense emission elsewhere in the broad-band filter. Such sources, with Lyman 
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Fig. 51. Night sky spectrum and the deployment of narrow band filters in 'quiet' 
regions corresponding to redshifted Lyman a emission as indicated below. The final 
optical window, corresponding to z ~6.9, was successfully exploited by lye et al 
(2006) to find two sources close to z ~7. 
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Fig. 52. The two-step process for locating high redshift Lyman a emitters (Hu et 
al 2004). (Left) Comparison of broad and narrow band magnitudes; sources with 
an unusual difference in the sense of being brighter in the narrow band filter rep- 
resent promising candidates. (Right) Spectroscopic follow-up reveals typically three 
possibilities - [O III] or [O II] at lower redshift, or Lya often characterized by its 
asymmetric line profile. 
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a equivalent widths in excess of several hundred A are interesting because 
they may challenge what can be produced from normal stellar populations. 




L(Lyo) [ergs s-'J 

Fig. 53. Comparison of the Lyman a luminosity functions at z=5.7 and 6.5 
from the surveys of Kashikawa et al (2006) and Shimasaku et al (2006); both 
spectroscopically-confirmed and photometric candidates are plotted. The decline in 
luminous emitters is qualitatively similar to trends seen over 3< z <6 in luminous 
continuum drop-outs. 

Malhotra & Rhoads (2004) were the first to consider the absence of evo- 
lution in the Lya LF as a constraint on the neutral fraction. Although they 
found no convincing change in the LF between z=5.7 and 6.5, the statistical 
uncertainties in both LFs were considerable. Specifically, below luminosities 
of L^ya ~10^^'^ ergs sec~^ no detections were then available. Hu et al (2005) 
have also appealed to the absence of any significant change in the mean Lya 



88 Richard S Ellis 



profile. Malhotra & Rhoads deduced the neutral fraction must be xhi <0.3 

at z ~Q supporting early reionization. However, Furlanetto et al (2005) rean- 
alyzed this constraint and indicated that strong emitters could persist even 
when Xhi —0.5 (c.f. Fig. 51). 

Kashikawa et al (2006) and Shimasaku et al (2006) have determined sta- 
tistically greatly improved Lyman a luminosity functions and discuss both 
spectroscopic confirmed and photometrically-selected emitters (Fig. 53). No 
decline is apparent in the abundance of low luminosity emitters as expected 
in an IGM with high xhi\ indeed the most significant change is a decline 
with redshift in the abundance of the most luminous systems. Although the 
change seems surprisingly rapid given the time interval is only 150 Myr, this 
is consistent with growth in the halo mass function (Dijsktra et al 2006). 

7.6 Results from Lensed Lycc Surveys 

The principal gain of narrow band imaging over other techniques in locating 
high redshift Lyman a emitters lies in the ability to exploit panoramic cameras 
with fields of view as large as 30-60 arcmin. Since cosmic lenses only magnify 
fields of a few arcmin or less, lensing searches are only of practical utility 
when used in spectroscopic mode. As discussed above, the gain in sensitivity 
can be factors of x30 or more, and given the small volumes explored, they 
are primarily useful in testing the faint end of the Lya luminosity function 
at various redshifts. A number of workers (e.g. Barkana & Loeb 2004) have 
emphasized the likelihood that the bulk of the reionizing photons arise from an 
abundant population of intrinsically-faint sources, and lensed searches provide 
the only practical route to observationally testing this hypothesis. 

A practical demonstration of a blind search for lensed Lya emitters is 
summarized in Figure 54. A long slit is oriented along a straight portion of 
the critical line (whose location depends on the source redshift). The survey 
comprises several exposures taken in different positions offset perpendicular 
to the critical line. Candidate emission lines are astrometrically located on 
a deep HST image and, if a counter-image consistent with the mass model 
can be located, a separate exposure is undertaken to capture both (as was 
the case in the source located by Ellis et al 2001). Unfortunately, continuum 
emission is rarely seen from a faint emitter and the location of a corresponding 
second image is often too uncertain to warrant a separate search. In this case 
contamination from foreground sources has to be inferred from the absence 
of corresponding lines at other wavelengths (Santos et al 2004, Stark et al 
2006b). 

Using this technique with an optical spectrograph sensitive to Lya from 
2.2< z <6.7, Santos et al (2004) conducted a survey of 9 lensing clusters and 
found 11 emitters probing luminosities as faint as LLya. — 10^° cgs, signifi- 
cantly fainter than even the more recent Subaru narrow band imaging searches 
(Shimasaku et al 2006). The resulting luminosity function is flatter at the faint 
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Fig. 54. Critical line mapping in AbeU 2218: how it works (Ellis et al 2001). The 
red curves show the location of the lines of very high magnification for a source at 
a redshift z—1 (dashed) and z—6 (solid). Blue lines show the region scanned at low 
resolution with a long-slit spectrograph. The upper right panel shows the detection 
of an isolated line astrometrically associated with (a) in the HST image for which a 
counter image (b) is predicted and recovered (see also inset to main panel). A higher 
dispersion spectrum aligned between the pair (yellow lines) reveal strong emission 
with an asymmetric profile in both (lower right panel). 



end than implied for the halo mass function and is consistent with suppression 
of star formation in the lowest mass halos. 

Stark et al (2006b) have extended this technique to higher redshift using 
an infrared spectrograph operating in the J band, where lensed Lya emitters 
in the range 8.5< z <10.2 would be found. This is a much more demanding 
experiment than that conducted in the optical because of the brighter and 
more variable sky brightness, the smaller slit length necessitating very precise 
positioning to maximize the magnifications and, obviously, the fainter sources 
given the increased redshift. Nonetheless, a 5a sensitivity limit fainter than 
10~^^ cgs corresponding to intrinsic (unlensed) star formation rates of 0.1 Af© 
yr~^ is achieved with the 10m Keck II telescope in exposure times of 1.5 hours 
per slit position. 

After surveying 10 clusters with several slit positions per cluster, 6 candi- 
date emission lines have been found and, via additional spectroscopy, it seems 
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Fig. 55. The volume density of sources of various star formation rates at z ~8-10 
required for cosmic reionization for a range of assumed parameters (blue hatched 
region) compared to the inferred density of lensed emitters from the survey of Stark 
at al (2006b). The open red symbol corresponds to the case if all detected emitters 
are at z ~10, the black symbols correspond to the situation if the two most promising 
candidates are at z ~10, and the dashed line corresponds to the 5a upper limit if 
none of the candidates is at z ~10. 



most cannot be explained as foreground sources. Stark et al estimate the sur- 
vey volume taking into account both the spatially-dependent magnification 
(from the cluster mass models) and the redshift-dependent survey sensitivity 
(governed by the night sky spectrum within the spectral band). 

Madau et al (1999) and Stiavelli et al (2004) have introduced simple pre- 
scriptions for estimating the abundance of star forming sources necessary for 
cosmic reionizations. While these prescriptions are certainly simple-minded, 
given the coarse datasets at hand, they provide an illustration of the implica- 
tions. 

Generally, the abundance of sources of a given star formation rate SFR 
necessary for cosmic reionization over a time interval At is 

Bum 



where B is the number of ionizing photons required to keep a single hy- 
drogen atom ionized, nn is the comoving number density of hydrogen at the 
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redshift of interest and fc is the escape fraction of ionizing photons. Figure 56 
shows the upshot of the Stark ct al (2006b) survey for various assumptions. 
The detection of even a few convincing sources with SFR ~0.1-1 Mq yr~^ in 
such small cosmic volumes would imply a significant contribution from feeble 
emitters at z ~10. Although speculative at this stage given both the uncer- 
tain nature of the lensed emitters and the calculation above, it nonetheless 
provides a strong incentive for continued searches. 

7.7 Lecture Summary 

In this lecture we have shown how Lyman a emission offers more than simply 
a way to locate distant galaxies. The distribution of line profiles, equivalent 
widths and its luminosity function can act as a sensitive gauge of the neutral 
fraction because of the effects of scattering by hydrogen clouds. Surveys have 
been undertaken using optical cameras and narrow band filters to redshifts 
z ~7. 

However, despite great progress in the narrow-band surveys, as with the 
earlier i-band drop outs, there is some dispute as to the evolutionary trends 
being found. Surprisingly strong evolution is seen in luminous emitters over 
a very short period of cosmic time corresponding to 5.7< z <6.5. And, to 
date there is no convincing evidence that line profiles are evolving or that 
the equivalent width distribution of emitters is skewed beyond what can be 
accountc;d for by normal young stellar populations. One suspects we will have 
to push these techniques to higher redshift which will be hard given the Lya 
line moves into the infrared where no such panoramic instruments are yet 
available. 

We have also given a brief tutorial on strong gravitational lensing. In about 
20 or so clusters, spectroscopic redshifts for sets of multiple images has enabled 
quite precise mass models to be determined which, in turn, enable accurate 
magnification maps to be derived. Remarkably faint sources can be found by 
searching along the so-called critical lines where the magnification is high. 
The techniques has revealed a few intrinsically faint sources and, possibly, the 
first glimpse of a high abundance of faint star forming sources at z ~10 has 
been secured. 
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8 Cosmic Infrared Background 

8.1 Motivation 

In this lecture we examine the role of cosmic backgrounds. First we make an 
important distinction. We are primarily concerned with extragalactic back- 
grounds composed of unresolved faint sources rather than the fairly isotropic 
microwave radiation that comes from the recombination of hydrogen. Such 
unresolved source background have played a key role in astrophysics. 

The pattern of discovery outside the optical and near-infrared spectral 
windows often goes like the following; the background is first discovered by 
sensitive detectors which do not have the angular resolution to see if there 
is any fine structure from faint sources. There is then some puzzlement as to 
its origin: for example, does it arise from an unforeseen population of sources 
beyond those already counted? In this manner, the X-ray background was 
identified as arising from active galactic nuclei and the sub-mm background 
from dusty star- forming galaxies. In each case we are concerned with sep- 
arating the counts of resolved sources to some detectability limit with the 
measured value of the integrated background. Key issues relate to the con- 
tribution of resolved sources and the removal of spurious (non-extragalactic) 
foreground signals. Excellent reviews of the subject have been provided by 
Leinert et al (1998), Hauser & Dwek (2001) and Kashlinksy (2005). 

In our case, we are interested in extending source counts of star forming 
galaxies beyond z ~6 and so the question of a near-infrared extragalactic 
background signal is of greatest relevance. The instruments concerned include 
cameras on Hubble Space Telescope and Spitzer SpaceTelescopc, but other 
space missions such as COBE and the Infrared Telescope in Space have pro- 
vided important results. 

8.2 Methodogy 

To understand the distinction between the resolved and unresolved compo- 
nents of the background, it is helpful to derive some fundamental relations for 
galaxy counts. 

In the magnitude system, the differential count slope 7 with increasing 
magnitude is 

dlogN{m) 
^ dm 

where N{m) is the differential number of galaxies per unit sky area (e.g. 
deg~^) in some counting bin dm. 

Now the contribution to the surface brightness of the extragalactic night 
sky from sources as a given magnitude m is: 

i^m) = l0-o-4(™+«°"''*) iV(m) 
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where the 0.4 factor arises from the relationship between flux and magni- 
tude. 

And the integrated surface brightness (the extragalactic background hght, 
EBL) is obtained by extending this integral to infinitely faint limits, viz: 

/oo 
ivim) dm 
-oo 

The bolometric equivalent of the EBL is then J ly dv. The EBL is often 
alternatively expressed as EBL = J v I^ dv. 

Surface brightness is usually expressed in nW m^^ steradian^^ although 
very early articles refer to some derivative of magnitudes deg~^. As infrared 
fluxes are often expressed in Janskies (1 Jansky = 10~^^ W m^^ Hz^^, a 
useful conversion is 1 nw m~^ ster~^ = 3000/A(/im) MJy ster^^. 

Examination of the above relation shows that if 7 >0.4, the surface bright- 
ness contribution from fainter sources outshines that of brighter ones and so 
the EBL will diverge. By deduction, therefore, the maximum contribution of 
resolved sources to the background will be where 7 ~0.4. If the slope of the 
counts turns down below 0.4 at some point, it may seem pointless to speculate 
that there is much information from fainter unresolved sources. In the case 
of searching for cosmic reionization however, where the flrst sources may be 
a distinct population, continuity in the source counts may not be expected. 
For this reason, the interesting quantity is the difference between a measured 
EBL and the integrated contribution from the faintest resolved sources. 




Fig. 56. (Left) Differential galaxy counts as a function of wavelength from the 
compilation of Madau & Pozzetti (2000). Note the absolute numbers have been 
arbitrarily scaled for convenience. (Right) Magnitude-dependent contribution of the 
counts to the surface brightness ii, of the extragalactic night sky. Depending on the 
waveband, the peak contribution occurs at uiab —20-25. 
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Madau & Pozzetti (2000) present a careful analysis of the deep optical 
and ncar-infrarcd counts at various wavelengths, mostly from Hubble Space 
Telescope data prior to the Ultra Deep Field (Figure 56). Extrapolation of the 
counts enables the contribution to the integrated light from known populations 
to be evaluated as a function of wavelength. The total EBL from this analysis 
is 55 nW m~^ sr~^ of which the dominant component lies longward of 1/zm. 
However, it must be remembered that galaxies are extended objects and so 
it is possible that significant light is lost in the outermost regions of each. As 
surface brightness is relativistically dimmed by the cosmic expansion, oc (1 + 
z)^, the contribution from distant sources could be seriously underestimated 
(Lanzetta et al 2002). 

In a similar fashion to the check we made that the integrated star formation 
history produces the present stellar mass density (Lecture 4), so it is possible 
to verify that the bolomctric output from stellar evolution should be consistent 
with the present mass density of stars. 

The bolometric radiation density pboi {t) is 



The check is a little bit trickier because a typical star formation history 
has to be assumed and presumably there is a large variety for different kinds 
of sources observed at various redshifts. The issue is discussed in detail by 
Madau & Pozzetti. 

8.3 Recent Background Measurements 

The goal for making extragalactic background light (EBL) measures is thus 
to determine the extent to which the measured value (or limit) exceeds that 
predicted from extrapolation of the galaxy counts. This might then provide 
some evidence for a new distant population such as the sources responsible 
for cosmic reionization. 

Various claims of an excess have been made in the 0.3 to 10/im wave- 
length region (Figure 57). Outside of this window the background appears 
to be entirely produced by known sources. Such EBL measurements are ex- 
tremely difficult to make for several reasons. An accurate absolute calibration 
is essential since the interesting signal is a 'DC difference' in surface bright- 
ness. The removal of spurious foregrounds is likewise troublesome: at some 
wavelengths, the foreground signal greatly exceeds the sought-after effect. 



Jo 

where ps is the star formation rate per comoving volume. 
The integrated EBL is then 
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Fig. 57. Summary of recent measures of the extragalactic background hght from 
the compilation of Dole et al (2006) . The cross-hatched region represents the region 
of claimed excess between the various background meaasures (labelled points) and 
the lower limits to the integrated counts (shown as upward arrows in the critical 
optical- near IR region). Interesting excesses are found in the 0.3 to 5 pim region. 

In the optical window, careful experiments have been undertaken by Mat- 
tila (1976), Dube et al (1979) and most recently Bernstein et al (2002). 
The most vexing foreground signals at optical wavelengths arise from airglow 
(emission in the upper atmosphere which has a time-dependent structure on 
fine angular scales), zodiacal light (scattering of sunlight by interplanetary 
dust which varies with the motion of the Earth along the ecliptic plane) and 
diffuse Galactic light or 'cirrus' (gas clouds illuminated by starlight). The 
wavelength dependence of these foregrounds (including the microwave back- 
ground itself which is considered a 'contaminant' in this respect!) is sum- 
marized in Figure 58. For the critical near-infrared region where redshifted 
light from early sources might be detected, airglow and zodiacal light are the 
dominant contaminants. 

Experiments have differed in the way these foregrounds might be removed. 
The use of a space observatory avoids airglow and observations at different 
Galactic latitudes can be used to monitor or minimize the effect of cirrus. 
Zodiacal light, the dominant foreground for all mid-infrared studies and all 
space-based near-infrared studies, is the most troublesome. 

To first order, the zodiacal light can be predicted depending on the geo- 
metrical configuration of the ecliptical plane, the earth, sun and target field. 
Although seasonal variations can be tracked and the spectrum of the zodiacal 
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Fig. 58. Wavelength dependence of the dominant foreground signals from the com- 
pilation by Kashlinsky (2005). Airglow and zodiacal light dominate at the wave- 
lengths of interest for stellar radiation from highly redshift sources. 

light is solar, imponderables enter such as the nature and distribution of the 
interplanetary dust and its scattering properties. 

The positive EBL excess measures at optical to near-infrared wavelengths 
have remained controversial. This is because all such experiments search for 
a 'DC signal' - a small absolute difference between two signals within which 
foregrounds have to be painstakingly removed. A brief study of the Bernstein 
et al (2002) experiment will make the difficulties clear. 

Bernstein et al (2002) claim a significant optical excess at 300, 550 and 
800nm from fields observed with Hubble's WFPC-2 and FOS. Airglow is elmi- 
nated since HST is above 90km and sources were removed in each HST image 
to Va_b=23.0. The zodiacal light spectrum was measured simultaneously with 
a ground-based optical telescope and iteratively subtracted from the HST 
data. 

At 550nm, the total WFPC-2 background measured after removing all 
sources was 105.7 ± 0.3 units (1 unit = 10~^ cgs ster"-'^ A ~-'^). The measured 
zodi background was 102.2 ±0.6 units. Galactic cirrus and faint galaxies be- 
yond yAB=23.0 were estimated as 0.8 and 0.5 units respectively. The claimed 
excess signal at this wavelength is 2.7±1.4 units. Not only is this only a 2 cr 
detection but it would only require a 2% error in the estimated Zodiacal light 
signal to be spurious. 

Studies of the infrared background advanced significantly with the launch 
of the COBE satellite which carried DIRBE - a 10 channel photometer oper- 
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ating in the 1-240 /ttm range with 0.7 degree resolution chopping at 32Hz onto 
an internal zero flux surface, and FIRAS - an absolute spectrometer with 7 
degree resolution operating in the 100 /xm to 5nim range. 

DIRBE measured the integrated background at 140 and 240/im using an 
elaborate timc-dcpcndcnt model of the Zodiacal light (Kelsall et al 1998, 
Hauser et al 1998). Schlegel et al (1998) combined these data with higher reso- 
lution IRAS 100/xm maps to improve removal of Galactic cirrus and to measure 
the temperature of the dust emission. Finkbcincr et al (2000) extended the 
model of Zodiacal light to provide the first detection at 60 and 100/um and 
Wright & Reese (2000) used 2MASS observations to improve Galactic source 
removal claiming a detection at 2.2/im. 

Dole et al (2006) have undertaken a very elegant analysis of 19000 Spitzer 
MIPS images. By centering the images on deep 24/im sources, they can eval- 
uate the statistical contribution of otherwise inaccessibly faint 70 and 160 fim 
sources. 

Matsumoto et al (2005) has extended these detections to shorter wave- 
lengths using a spectrometer onboard the Infrared Telescope in Space (IRTS). 
This signal has the tantalising signature of a distant star-forming stellar pop- 
ulation - a steeply-rising continuum down to 1/um and a discontinuity when 
extended to include optical data (Figure 57). Mattila (2006) has argued this 
signal represents an artefact arising from incorrect foreground removal. He 
has likewise criticised the Bernstein et al (2002) optical detections (Mattila 
2003). 

Supposing the DIRBE/2MASS/IRTS detections to be real, what would 
this imply? Assuming the J-band background (>2.5 nW m~^ ster~^) arises 
from z ~9 Lya emission, Madau & Silk (2005) calculate the associated stellar 
mass that is produced. They find an embarrassingly high production rate, 
corresponding to J7* = 0.045 Qb', in other words almost all the stars we see 
today would need to be produced by z ~9 to explain the signal. Likewise 
the ionizing flux produced by such star formation would be in excess of that 
required to explain the WMAP optical depth. 

8.4 Fluctuation Analyses 

Kashlinksky et al (2002, 2005) have argued that the difficulties inherent in 
extracting the EBL signal by the DC method may be alleviated by considering 
the angular fluctuation spectrum 5F{6) and its 2-D Fourier transform P^in)- 
In more detail, the fluctuation in the measured background is 

5F{e) = F{9)- <F>= (271-)"^ y" 6Fqexp{=iq.9)cfq 

and 



P2{q) =< \dF,f> 
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The success of the method depends on whether the various foreground con- 
tributions to P2{q) can be readily distinguished from one another. Although 
advantageous in using independent information from the DC measurements, 
it is unfortunately not easy to intercompare the experiments or to interpret 
the fluctuation analyses. 

Kashlinksky et al (2005) recently applied this method to deep IRAC im- 
ages. Sources were extracted to reasonably faint limits (0.3 /iJy or niAB —22- 
25) . The pixels associated with each image were masked out in a manner that 
could be adjusted depending on the significance of the source over the back- 
ground and its area. The data was split into two equal halves (A, B) and the 
power spectrum of the signal {Ps{q), S=A+B) was compared to that of the 
noise (P/v((7), N=A-B). Figure 59 illustrates the residual signal, Ps{q)—PN{q), 
and the fluctuations {q'^P2{q)/2i:)2 in 3 IRAC fields as a function of angular 
scale (27r/g arcsec) in the four IRAC channels. Although the positive fluctu- 
ations on small scales are consistent with shot noise from the galaxy counts 
(solid lines in both panels), the excess is particularly prominent on scales of 
1-2 arcmin and consistently in fields of various Galactic latitudes and in all 
four IRAC bandpasses. 




Fig. 59. Detection of fluctuations in the infrared background from the analysis by 
Kashlinky et al (2005). The upper panel shows the angular power spectrum as a 
function of scale for 3 IRAC fields plotted for all four IRAC channels (3.6 through 
8 ^m). The solid line shows the effect of shot noise from the galaxy counts. This 
fits the fiuctuations on small scales but there is a significant and consistent excess 
on large scales. Bottom panels refer to fluctuations for one IRAC fleld plotted as 
(?^-P2(q)/27r) 2 with the predicted contributions from various populations of unre- 
solved distant sources shown as dotted lines. 
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Kashlinksky et al argue that the excess signal is fairly flat spectrally in v ly 
ruling out any instrumental or Galactic contamination. The constancy of the 
signal over the wide range in Galactic latitude likewise suggests the signal is 
extragalactic. However, zodiacal light is difficult to eliminate as a contaminant. 
Although the signal persists when the Earth is 6 months further around in its 
orbit, no analysis of the likely zodi fluctuation spectrum is presented. 

If the signal persists, what might it mean for sources fainter than 0.3 /iJy? 
Without a detailed redshift distribution, it is hard to be sure. The cluster- 
ing pattern is remarkably strong given the depth of the imaging data (much 
deeper fields have subsequently been imaged but no independent analyses 
have yet been published). In summary, this is an intriguing result, somewhat 
unexplained and in need of confirmatior(f| 



8.5 EBL Constraints from TeV Gamma Rays 

High energy (TeV) gamma rays are absorbed in the earth's atmosphere and 
converted into secondary particles forming an 'air shower'. Cerenkov light is 
generated - a beam of faint blue light lasting a few xlO~^ sec is produced 
illuminating an area of 250m in diameter on the ground. 

The significance of searching for TeV gamma rays is that they interact 
with the sea of 1-10 /im (infrared background) photons via the pair-creation 
process, viz: 

77 — > e"*" e~ 

producing attenuation in distant sources (blazars) whose spectral energy 
distributions are assumed to be power laws. 

dN r 

oc 

dE 

The strength of the attention, measured as a change in F, for the most 
distant accessible blazars is thus a measure of the degree of interaction between 
the gamma rays and the ambient infrared background. 

The HESS team (Aharonian et al 2006) have recently analyzed the gamma 
ray spectrum of a blazar at 2=0.186 and fitted its energy spectrum. As a result 
they can predict an upper limit to the likely infrared background for various 
assumptions (Figure 60). The constraint is most useful in the wavelength 
range 1-10/im where the associated gamma ray spectral data is of high quality 
(top axis of Fig. 60). This new method is intriguing. Although in detail the 
constraint is dependent on an assumed form for the TeV energy spectrum 
of blazars, for all reasonable assumptions the acceptable background is much 
lower at 1-4 yitm than claimed by the traditional experiments. 



Recent developments are discussed by Sullivan et al 2006 and Kashlinsky et al 
2006 
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Fig. 60. Constraints on the cosmic infrared background from the degree of atten- 
uation seen in the high energy gamma ray spectrum of a distant blazar (Aharonian 
et al 2006). The black curves indicate likely upper limits to the infrared background 
for various assumptions. This independent method argues the near-IR detections by 
Matsumoto et al (2005, blue circles) and others are spurious. 

8.6 Lecture Summary 

In this lecture we have learned that extragalactic background light measure- 
ments in principle offer an important and unique constraint on undiscovered 
distant source populations. 

However, the observations remain challenging and controversial because 
of instrumental effects, related to the precision of absolute calibrations, and 
dominant foregrounds. In the 1-5/im spectral regions positive detections in 
excess of extrapolations of the source counts have been reported. This is in- 
teresting because this is the wavelength range where an abundant and early 
population of faint, unresolved z >10 sources would produce some form of 
infrared background. 

At the moment, the claimed detections at 1-5 /zm appear unreasonably 
strong. Too many stars would be produced by the star formation associated 
with this background signal and recent ultrahigh energy spectra of distant 
blazars likewise suggest a much weaker infrared background. 
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IRAC is a particularly promising instrument for the redshift range 10< 
z <20 and fluctuation analyses offer a valuable independent probe of the 
background, although it is hard to interpret the results and compare them 
with the more traditional DC level experiments. 

Despite the scemlingly chaotic way in which the subject of the cosmic 
extragalactic background light has unfolded in the past decade, it is worth 
emphasizing that the motivation remains a strong one. Any information on the 
surface brightness, angular clustering and redshift range of radiation beyond 
the detected source counts will be very valuable information ahead of the 
discoveries made possible by future generations of facilities. 
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9 Epilogue: Future Prospects 

9.1 Introduction 

The students attending this lecture course are hving at a special time! During 
your lifetime you can reasonably confidently expect to witness the location of 
the sources responsible for cosmic reionization, to determine the redshift range 
over which they were active and perhaps even witness directly the 'cosmic 
dawn' when the first stellar systems shone and terminated the 'Dark Ages' ! 

We have made such remarkable progress in the past decade, reviewed here, 
that such a bold prediction seems reasonable, even for a cautious individual 
like myself! We have extended our fundamental knowledge of how various 
populations of galaxies from 0< z <3 combine to give us a broad picture 
of galaxy evolution, while extending the frontiers to z ~7 and possibly be- 
yond. Certainly many issues remain, including the apparent early assembly of 
certain classes of quiescent galaxy, the abundance patterns in the intergalac- 
tic medium and the apparent 'downsizing' signatures seen over a variety of 
redshifts. However, the progress has been rapid and driven by observations. 
Accompanying this is a much greater synergy between theoretical predictions 
and observations than ever before. 

We have seen that the question of 'First Light' - the subject of this course 
- is the remaining frontier for observations of galaxy formation. The physical 
processes involved are poorly understood and thus observations will continue 
to be key to making progress. In this final lecture, I take out my crystal ball 
and consider the likely progress we can expect at optical and near-infrared 
wavelengths with current and near-term facilities ahead of those possible with 
the more ambitious new observatories such as the James Webb Space Tele- 
scope (JWST) and a new generation of extremely large ground-based tele- 
scopes. 

9.2 The Next Five Years 

Panoramic imaging with large optical and near-infrared cameras on telescopes 
such as Subaru, UKIRT and VISTA will enable continued exploration of the 
abundance of luminous drop-outs and Lya emitters over 5< z <7, reducing 
the currently troublesome issues of cosmic variance. Deeper data will continue 
to be provided from further fields taken with ACS and NICMOS. With further 
spectroscopic surveys on large telescopes, we can expect improved stellar mass 
density estimates at z ~5-6 (Eyles et al 2006). However, strong gravitational 
lensing may still be the only route to probing intrinsically fainter sources, 
particularly beyond z ~7. 

More detailed characterization of the properties of the most massive galax- 
ies at z ~5-6 will also be worthwhile, in addition to continuing to deduce sta- 
tistical properties such as abundances and luminosity functions.. Laser-guide 
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Fig. 61. Detailed studies of a lensed z ~5 galaxy (Swinbank et al 2006). (Top) Re- 
constructed color HST VI image of the z=4.88 arc in the cluster RCS0224-002. The 
inset shows the effect of 0.8 arcsec seeing on the reconstruction, thereby demonstrat- 
ing the advantage of lensing. In the source plane, the galaxy is 2.0x0.8 kpc. (Bottom) 
[O II] velocity field obtained during a 12 hr VLT SINFONI exposure without AO. 
Spatial comparison with the Lya field gives clear evidence of significant bi-polar 
outflows. 



star adaptive optics is now in widespread use on our 8-10 meter class tele- 
scopes and, with integral field spectrographs, can be used to probe the resolved 
dynamics of distant galaxies (Genzel et al 2006). Application of these tech- 
niques to the most intense star-forming lensed systems at z >4 (Swinbank et 
al 2006, Figure 61) is already providing unique insight into the physical state 
of the earliest galaxies. Extending the same techniques with adaptive optics 
will advance progress and address process on remarkably small scales (~200 
pc). 
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A number of special purpose ground-based instruments are also being de- 
veloped, both with and without adaptive optics, to probe for z >7 dropouts 
and Lyman a emitters. These include: 

• The Dark Age Z Lyman Explorer (DAZLE, Horton et al 2004 j^, a mod- 
erate field (7 arcmin) infrared narrow-band imager with airglow discrim- 
ination whose goal is to reach a sensitivity of ~ 10~^^ cgs in a night of 
VLT time. This corresponds to a limiting star formation rate of ~ IA/q 
yr^^ in quiet regions of the airglow spectrum at z c^ij .1 and 9.9. Given the 
modest field and the need for dedicated spectroscopic follow-up, a long 
term campaign is envisaged to conduct a reliable census. 

• The Gemini Genesis Survey (GGS 0) which uses F2T2, a clone of the 
tunable filter destined for JWST (Scott et al 2006). This instrument is 
planned to work behind the Gemini Observatory's multi-conjugate adap- 
tive optics facility and image, in steps of wavelength, regions lensed by 
foreground clusters. The success of such a strategy will depend on both 
the angular size and line widths of z >7 Lyman a emitters and the abun- 
dance of intrinsically faint examples. Ellis et al (2001) found the lensed 
emitter at z ~5.7 was less than 30 milli-arcsec across with a line width 
of only 100-200 kms sec""'^. If such tiny emitters are the norm at z >7, 
the GGS may go much deeper than extant spectroscopic lensed searches 
(Stark et al 2006b). 

The exciting redshift range 7< z < 12 will also be the province of improved 
drop-out searches using the instrument WFC3 slated to be installed on Hubble 
Space Telescope in 2008-JO- The infrared channel of this instrument spans 
850 - 1170 nm with a field of ~ 2 arcmin at an angular resolution of 0.13 
arcsec pixel" ^. This resolution is coarser than that of adaptive optics- assisted 
instruments such as F2T2. The principal gain over ground-based instruments 
will be in deep broad-band imaging free from airglow. The survey efficiency 
is about an order of magnitude better than that of NICMOS. WFC3 also has 
two infrared grisms which will be helpful in source discrimination. 

A major stumbling block at the moment, even at z ~5-6, is efficient spec- 
troscopic follow-up of dropout candidates. As we discussed in §6, photometric 
redshifts have unfortunately become de rigeur in statistical analyses of lu- 
minosity densities and luminosity functions (Bouwens et al 2006), yet their 
precision remains controversial. For z- and J-band dropous beyond z ~7, pho- 
tometric redshifts will be even less reliable. Spitzer detections will be harder 
and fewer, and the typical source may have only 2-3 detected bands. Candi- 
dates may be found in abundance but how will they be confirmed? 



tLttp://www.ast.cam.ac.uk~optics/dazle/| 

http://odysseus. astro. utoronto.ca/ggs-blog/?page_id=2 
http://www.stsci.edu/hst/wfc3l 
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The various 8-lOm telescopes are now building a new generation of cryo- 
genic near-infrared multi-slit spectrographs. It can be hoped that long expo- 
sures with these new instruments will be sufficient to break this impasse. 

9.3 Beyond Five Years 

A number of facilities are being planned, motivated by the progress discussed 
in these lectures. These include: 




Fig. 62. The 6.5m James Webb Space Telescope: Then (2013 in orbit) and and 
Now (2005, fuU scale model). 



• James Webb Space Telescope - a 6.5 meter optical- infrared space observa- 
tory (Gardner et al 2006) for which a primary mission is the detailed study 
of the earliest star-forming galaxies (Figure 62). The facility is currently 
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Fig. 63. The proposed US-Canadian Thirty Meter Telescope (www.tnit.org) now 
in the detailed design phase. 

planned to have a near-infrared imager NIRCam, a spectrograph NIRSpec 
with both integral field and limited multi-object modes, a tunable filter (as 
F2T2) and a mid-IR imager, MIRI. StiavcUi (2002) has presented a cogent 
summary of the likely strategies of using this facility, due to be launched 
in 2013, for studies of the earliest systems. 
• Extremely Large Telescopes on the ground including the US-Canadian 
Thirty Meter Telescope (TMT, see www.tmt.org, Figure 63) which will 
have a diffraction- limited near- infrared imager/spectrograph (IRIS) and a 
adaptive-optics assisted spectrograph (IRMOS) with multiple deployable 
integral field units - area mapping units which can be arranged not only 
to scan regions surrounding luminous sources but also to undertake multi- 
object studies in crowded regions (Figure 64). At the time of writing, TMT 
is expected to be operational from 2016 onwards. 

These will complement redshifted 21cm line tomography with radio fa- 
cilities and address the key questions of the escape fraction of photons from 
star-forming sources and how they create ionized bubbles which merge to 
cause reionization. 

It is quite likely that, by 2013, the redshift range containing the earliest 
galactic sources, estimated at present to be 10< z <20 perhaps, will have been 
refined sufficiently by special-purpose instruments on our existing 8-10 meter 
class telescopes. Thus one can surmise that that both JWST and future ELTs 
will be used for much more challenging work related to the physical process 
of reionization, as well as the chemical maturity of the most luminous sources 
found at high redshift. 
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Fig. 64. Flexibility in the deployment strategy for the multiple integral field units 
(yellow squares) for the proposed TMT infrared multi-object adaptive optics assisted 
instrument IRMOS. The IFUs can be distributed around the full field of 2 arcmin 
in classic multi-object mode. Alternatively, the IFUs can be configured together to 
map a small 6x6 arcsec field at high angular resolution in 'blind' mode. Such 
fiexibility will be important in mapping Lya emission at high redshift in various 
situations. 

An obvious partnership between JWST and TMT, for example, which 
would complement the 21cm studies, would be to (i) search for the extent and 
topology of faint Lya emission in ionized bubbles around JWST-selected lu- 
minous star forming galaxies and, (ii) pinpoint early sources for spectroscopic 
scrutiny so as to identify signatures of Population III stars. 

At the present time there are so many imponderables in our knowledge of 
the earliest sources, that even the design parameters for the ELT instruments 
is a considerable challenge. How big are the faintest Lya emitters? What are 
the typical line widths in km sec~^? How big are the ionized bubbles at a 
given redshift? And, crucially, what is the surface density of various types 
of star- forming galaxies. Flexibility in the design of survey strategies will be 
crucial with instruments such as TMT's IRMOS (Fig. 64). 

Any information we can glean on the properties of z 2±10 sources in the 
next 5 years will be valuable in optimizing how to move forward when these 
magnificent new generation telescopes are made available to us. Indeed, it is 
foolhardy to wait! Time and again, we can retrospectively look back at what 
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we thought we would accomplish with our planned facilities and we always 
find that we achieved more than we expected! 
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